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Summary 

Pratt & Whitney has developed a Broadband Fan Noise Prediction System (BFaNS) for turbofan 
engines. This system computes the noise generated by turbulence impinging on the leading edges of the 
fan and fan exit guide vane, and noise generated by boundary-layer turbulence passing over the fan 
trailing edge. BFaNS has been validated on three fan rigs that were tested during the NASA Advanced 
Subsonic Technology Program (AST). The predicted noise spectra agreed well with measured data. The 
predicted effects of fan speed, vane count, and vane sweep also agreed well with measurements. 

The noise prediction system consists of two computer programs: SetupBFaNS and BFaNS. 
SetupBFaNS converts user-specified geometry and flow-field information into a BFaNS input file. From 
this input file, BFaNS computes the inlet and aft broadband sound power spectra generated by the fan and 
FEGV. The output file from BFaNS contains the inlet, aft and total sound power spectra from each noise 
source. 

This report is the first volume of a three-volume set documenting the Broadband Fan Noise Prediction 
System: 


Volume 1: Setup BFaNS User’s Manual and Developer’s Guide 
Volume 2: BFaNS User’s Manual and Developer’s Guide 
Volume 3: Validation and Test Cases 

The present volume begins with an overview of the Broadband Fan Noise Prediction System, 
followed by step-by-step instructions for installing and running Setup BFaNS. It concludes with technical 
documentation of the Setup BFaNS computer program. 


iii 
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Broadband Fan Noise Prediction System for Turbofan Engines 

Volume 1: Setup BFaNS User’s Manual and Developer’s Guide 


Bruce L. Morin 
Pratt & Whitney 

East Hartford, Connecticut 06108 


1.0 Introduction 

In early turbofan engines, the fan blade passage tone and higher harmonics dominated the fan noise 
level when compared to fan broadband noise. However, engine makers have been very successful at 
reducing fan tone noise by decreasing fan tip speed, providing ample spacing between airfoil rows, 
carefully selecting the airfoil counts, and acoustically treating the fan duct. In fact, tone noise reduction 
has been so great that fan broadband noise is now a major contributor to the overall engine noise level. 

Under NASA funding, Pratt & Whitney has developed a Broadband Fan Noise Prediction System 
(BFaNS) for turbofan engines. This prediction system is built around acoustic theories developed by 
Hanson (Refs. 1 to 4) and Glegg (Refs. 5 to 7). These theories account for noise generated by turbulence 
impinging on the leading edges of the fan and fan exit guide vane (FEGV), and noise generated by 
boundary-layer turbulence passing over the fan trailing edge. 

Figure 1 shows the fan broadband noise sources that can be calculated with BFaNS. The blade tip- 
noise and blade self-noise calculations are limited to subsonic blade-tip relative Mach numbers. There are 
no limitations on the other sources. 

This report begins with an overview of the Broadband Fan Noise Prediction System, followed by 
step-by-step instructions for installing and running Setup BFaNS. It concludes with technical 
documentation of the SetupBFaNS computer program. Volume 2 of this report series provides 
instructions for running BFaNS, and Volume 3 provides validation and test cases. 

2.0 General Information 

This section provides general information about the Broadband Fan Noise Prediction System, 
including the orientation of the coordinate system, the velocity and angle conventions, and the 
conventions used to describe the flow path and airfoil geometries. 


Interaction of 
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Figure 2. — Organization of BFaNS. 




Rear View 
Figure 3. — Engine coordinate system. 


2.1 Overview 

Figure 2 shows the organization of the Broadband Fan Noise Prediction System, which consists of two 
computer programs: SetupBFaNS and BFaNS. SetupBFaNS converts user-provided geometry and flow- 
field information into a BFaNS input file (bfans. input). From this input file, BFaNS computes the inlet and 
aft broadband sound power spectra generated by the fan and FEGV. The output file (bfans. output) contains 
the inlet, aft and total sound power spectra from each noise source shown in Figure 1, along with the 
combined sound power spectra from all sources. Both Setup BFaNS and BFaNS use control files (also 
referred to as author files) that allow the user to modify the way the programs operate. 

2.2 Engine Coordinate System 

Figure 3 shows the engine coordinate system used in Setup BFaNS and BFaNS. The z-direction is 
measured along the engine axis with z-values increasing from the inlet toward the exit. The axial location of 
the origin is arbitrary (for convenience, it is shown at the blade stacking line). When viewed from the rear, 
the x-axis is horizontal with values increasing toward the right. Likewise, the y-axis is vertical with values 
increasing in the upward direction. The corresponding cylindrical coordinate system is defined as follows: 

r - -Jx 2 + y 2 

tan 9 = — 

x 

Defined in this manner, 9 increases in the counter-clockwise direction when viewed from the rear. 

2.3 Velocity Conventions 

Figure 4 shows the velocity conventions used in Setup BFaNS and BFaNS. The conversion from the 
stationary reference frame to the rotating reference frame is: 

W = C-U 
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RearView 


Figure 4. — Velocity conventions. 


From this conversion formula, we get the following relationships: 

W Z =C Z 

w r = c r 

W e = C Q -U 

Note that U is positive when the fan rotates counter-clockwise viewed from the rear, and negative when 
the fan rotates clockwise viewed from the rear. The meridional velocity component (which is the same in 
both reference frames) is defined as: 


c =w = Jr 2 + c 2 

'' m \ y - / z ' 


2.4 Flow-Angle Definitions 

The flow angles are computed from: 


tan|3 = 


Wq_ 

w m 


tana = 


Ql 

r 

y-'m 


tancj) = 

C 

'-'Z 


El 

^z 


The a and (3 angles are measured from the meridional direction to facilitate stripwise application of the 
acoustic theories used in BFaNS (see Volume 2: BFaNS User’s Manual and Developer’s Guide). 


2.5 Flow-Path Coordinates 

Figure 5 shows the side view of a fan stage, along with some of the important terminology used in 
Setup BFaNS and BFaNS. The flow path is axisymmetric and consists of the hub surface, shroud surface, 
upper splitter surface and lower splitter surface. These surfaces must be defined in terms of their ( r,z ) 
coordinates. Specific points of interest are: the axial location of the hilite, the nose-cone leading edge, the 
splitter leading edge, the duct exit, and the blade and vane trailing-edge tips. Stations 1, 2, 3, and 4 
indicate the locations where flow-field information is available (via experiment or prediction). Flow-field 
information must be available upstream from the blade leading edge, downstream from the blade trailing 
edge (but upstream from the splitter leading edge), upstream from the vane leading edge (but downstream 
from the splitter leading edge), and downstream from the vane trailing edge. 
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nose cone 
leading edge 

Figure 5. — Flow-path terminology. 


rQ* 



Figure 6. — Airfoil coordinate system. 


2.6 Airfoil Coordinates 

Figure 6 shows the coordinate system used to define airfoil geometry. The leading and trailing edge 
metal angles are measured clockwise from the suction-side tangential; therefore, the metal angles are 
always positive. The asterisks on x, y and 0 indicate that the x*y* -plane ma y be rotated relative to the 
xy-planc shown in Figure 3. Likewise, the asterisk on z indicates that the axial location of the x*y* -plane 
may be translated relative to the xy-plane shown in Figure 3. Defining the airfoil in this manner allows the 
blade and vane geometry to be defined independently from each other, and independently from the flow 
path. Section 4.2.2 describes the coordinate transformations used by Setup BFaNS to ensure that the 
blade and vane are properly located in the flow path, and are properly oriented relative to the direction of 
rotation. 
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2.7 Sweep and Lean Definitions 


Figure 7 shows a side view and rear view of an airfoil. In SetupBFaNS, the leading edge is described 
parametrically as a function of radius (z e and 0 te versus r). Therefore, the sweep and lean are computed 
from, 


tanij/^ = 


dzie 

dr 


tan\|/^ = r 


dQ/ e 

dr 


Both of these quantities depend on radius. The overall sweep and lean are computed from: 

A z u 


tanTfr = 


A r 


tan T/ = r ^ le 
A r 

where the A is applied from the hub to the tip, and r is the mean radius. Sections 4.8.5 and 4.8.6 provide 
details of the parametric description for the airfoil leading edge. 



(Side View) 



(Side View) 



(Rear View) 



(Rear View) 


Figure 7. — Sweep and lean conventions. 


N ASA/CR— 20 1 0-2 1 6898/VOL 1 


5 



3.0 Installation and Execution of SetupBFaNS 

This section provides instructions for installing and running Setup BFaNS. Refer to Section 4.0 for 
technical documentation of the computer program. 

3.1 Installation 

0 provides instructions for installing the Broadband Fan Noise Prediction System (Setup BFaNS and 
BFaNS) on Unix-based computers. Both computer programs are written in FORTRAN 77. Appendix C 
explains how to run the Setup BFaNS test case. 

3.2 Run Procedure 

Setup BFaNS converts user-specified geometry and flow-field information into a BFaNS input file 
(bfans. input). Table 1 shows the step-by-step procedure for running Setup BFaNS. Sections 3.2.1 through 
3.2.10 describe the details of each step. 


TABLE 1.— RUN PROCEDURE FOR SETUP BFANS 


Step 1 

Create directory structure 

Step 2 

Create flow-path geometry file 

Step 3 

Create blade and vane geometry files 

Step 4 

Create flow-field files 

Step 5 

Create Setup BFaNS input file 

Step 6 

Create Setup BFaNS author file (optional) 

Step 7 

Execute Setup BFaNS 

Step 8 

Examine boundary-layer characteristics 

Step 9 

Examine wake characteristics 

Step 10 

Examine output file 


3.2.1 Create Directory Structure 

The user must create a directory that contains the following sub-directories: 

• A sub-directory that contains the files that describe the geometry of the flow path, the blade, and 
the vane. The default name of this sub-directory is bfans_geometry, but the user can over-ride 
this name in the Setup BFaNS author file (see Section 3.2.6, Card 5002). 

• (Optional) a sub-directory that contains the files that describe the inviscid streamline -predicted 
flow field near the blade leading edge, blade trailing edge, vane leading edge, and vane trailing 
edge for each configuration and operating point. The default name of this sub-directory is 
bfans_sline, but the user can over-ride this name in the Setup BFaNS author file (see Section 
3.2.6, Card 5004). 

• (Optional) a sub-directory that contains the files that describe the viscous CFD-predicted flow 
field near the blade leading edge, blade trailing edge, vane leading edge, and vane trailing edge 
for each configuration and operating point. The default name of this sub-directory is bfans_cfd, 
but the user can over-ride this name in the Setup BFaNS author file (see Section 3.2.6, Card 
5006). 

• (Optional) a sub-directory that contains the files that describe the measured flow field near the 
blade leading edge, blade trailing edge, vane leading edge, and vane trailing edge for each 
configuration and operating point. The default name of this sub-directory is bfans_data, but the 
user can over-ride this name in the Setup BFaNS author file (see Section 3.2.6, Card 5008). 

• A sub-directory (or several sub-directories) where Setup BFaNS will be run. These sub- 
directories are referred to as working sub-directories. Typically, separate working sub-directories 
are created for each operating condition (e.g., approach, cutback, and sideline). 
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3.2.2 Create Flow-Path Geometry File 

The user must create a file that contains the ( r,z ) coordinates of the hub, the shroud, the upper splitter, 
and the lower splitter surfaces (see Figure 5). The axial coordinates can be measured from any reference 
point, though the same reference point must be used for each surface. The length of the filename must be 
less than or equal to 30 characters. This file name is specified in the Setup BFaNS input file (see Section 
3.2.5, block 1). The flow-path geometry file must be located in bfans_geometry or the directory specified 
in the Setup BFaNS author file (see Section 3.2.6, Card 5002). 

Table 2 explains the file format for the flow-path geometry file. Appendix D shows a sample file 
along with a section of FORTRAN code used to read the file. Note that this file contains 27 “blocks” of 
input. Throughout this report, a “block” refers to one or more variables entered on a single line of input 
(e.g., block 1), or a single array variable entered on several lines of input (e.g., block 4). 


TABLE 2.— FORMAT OF FLOW-PATH GEOMETRY FILE 


Block 

Variable Name 

Type 

Format 

Description 

1 

filedescr 

character 

a50 

description of data in blocks 2-6 

2 

nduct h 

integer 

free 

number of points along hub 

3 

parmname 

character 

a50 

parameter name for block 4 

4 

zhub (i) 
i=l, nduct h 

real 

free 

hub axial coordinates (in) 

5 

parmname 

character 

a50 

parameter name for block 6 

6 

rhub ( i ) 
i=l, nduct h 

real 

free 

hub radial coordinates (in) 

7 

blank 

character 

al 

blank line 

8 

filedescr 

character 

a50 

description of data in blocks 9-13 

9 

nsplt u 

integer 

free 

number of points along upper splitter 

10 

parmname 

character 

a50 

parameter name for block 11 

11 

zsup (i) 
i=l, nsplt u 

real 

free 

upper splitter axial coordinates (in) 

12 

parmname 

character 

a50 

parameter name for block 13 

13 

rsup (i) 
i=l, nsplt u 

real 

free 

upper splitter radial coordinates (in) 

14 

blank 

character 

al 

blank line 

15 

filedescr 

character 

a50 

description of data in blocks 16 - 20 

16 

nsplt 1 

integer 

free 

number of points along lower splitter 

17 

parmname 

character 

a50 

parameter name for block 18 

18 

zslo (i) 
i=l, nsplt 1 

real 

free 

lower splitter axial coordinates (in) 

19 

parmname 

character 

a50 

parameter name for block 20 

20 

rslo (i) 
i=l, nsplt 1 

real 

free 

lower splitter radial coordinates (in) 

21 

blank 

character 

al 

blank line 

22 

filedescr 

character 

a50 

description of data in blocks 23 - 27 

23 

nduct t 

integer 

free 

number of points along shroud 

24 

parmname 

character 

a50 

parameter name for block 24 

25 

ztip (i) 
i=l, nduct t 

real 

free 

shroud axial coordinates (in) 

26 

parmname 

character 

a50 

parameter name for block 27 

27 

rtip (i) 
i=l, nduct t 

real 

free 

hub radial coordinates (in) 


3.2.3 Create Blade and Vane Geometry Files 

The user must create a file that contains the (x*,y*,z*) coordinates of the leading and trailing edges of 
the blade, along with the corresponding metal angles (see Section 2.6 for details). The user must create a 
similar file for the vane. The axial coordinates can be measured from an arbitrary reference point (the 
reference can be different for the blade and the vane, and can also be different than the reference used to 
describe the flow path). The length of the filenames must be less than or equal to 30 characters. These file 
names are specified in the Setup BFaNS input file (see Section 3.2.5, blocks 8 and 11). These files must 
be located in bfans_geometry or the directory specified in the Setup BFaNS author file (see Section 
3.2.6, Card 5002). 
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Table 3 explains the format of the blade and vane geometry fdes. Appendix E and Appendix F show 
sample blade and vane geometry files along with sections of FORTRAN code used to read these files. 


TABLE 3.— FORMAT OF BLADE AND VANE GEOMETRY FILES 


Block 

Variable Name 

Type 

Format 

Description 

1 

filedesc 

character 

a50 

file description 

2 

numrle 

integer 

free 

number of points defining blade or vane 

3 

parmname 

character 

a50 

parameter name for block 4 

4 

xle (i) 
i=l , numrle 

real 

free 

l.e. x-coordinates (in) 

5 

parmname 

character 

a50 

parameter name for block 6 

6 

yle (i) 
i=l , numrle 

real 

free 

l.e. y-coordinates (in) 

7 

desc 

character 

a50 

parameter name for block 8 

8 

zle (i) 
i=l , numrle 

real 

free 

l.e. z-coordinates (in) 

9 

parmname 

character 

a50 

parameter name for block 10 

10 

betale (i) 
i=l , numrle 

real 

free 

l.e. metal angle (deg) 

11 

parmname 

character 

a50 

parameter name for block 12 

12 

xte (i) 
i=l , numrle 

real 

free 

t.e. x-coordinates (in) 

13 

parmname 

character 

a50 

parameter name for block 14 

14 

yte (i) 
i=l , numrle 

real 

free 

t.e. y-coordinates (in) 

15 

parmname 

character 

a50 

parameter name for block 16 

16 

zte (i) 
i=l , numrle 

real 

free 

t.e. z-coordinates (in) 

17 

parmname 

character 

a50 

parameter name for block 18 

18 

betate (i) 
i=l , numrle 

real 

free 

t.e. metal angle (deg) 


3.2.4 Create Flow-Field Files 

The user must create files that describe the flow field at the following four locations (see Figure 5): 

• Station no. 1 — Upstream from the blade leading edge 

• Station no. 2 — Downstream from the blade trailing edge, but upstream from the splitter leading 
edge 

• Station no. 3 — Upstream from the vane leading edge, but downstream from the splitter leading 
edge 

• Station no. 4 — Downstream from the vane trailing edge 

The flow-field information can be provided by any combination of predictions (e.g., an inviscid 
streamline deck, or viscous CFD) or experimental measurements. Table 4 explains the file format used for 
viscous CFD and measured data, and Appendix G shows the section of FORTRAN code used to read 
these files. 

and 0 show similar information for inviscid streamline predictions. 

The length of each filename must be less than or equal to 30 characters. These names are specified the 
SetupBFaNS input file (see Section 3.2.5, blocks 13, 15, 17 and 19). The flow-field files must be located 
in bfans_sline, bfans_data, bfans_cfd or the directories specified in the Setup BFaNS author file (see 
Section 3.2.6, Cards 5004, 5006, or 5008). 

For viscous CFD predictions and experimental data, the axial coordinates must be measured from the 
blade trailing-edge tip. The blade trailing-edge tip was chosen as the reference point because of its 
convenience for experimentally measured data. For inviscid streamline predictions, the axial coordinates 
must be measured from the same reference that is used when defining the flow path. 
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The velocity, total pressure and total temperature can be specified at different locations. This feature 
is built into the program because the total-pressure and total-temperature measurements are often taken at 
a different axial location than the velocity measurements (e.g., see Volume 3: Validation and Test Cases). 
The rotational speed must be consistent with the direction of rotation specified in the SetupBFaNS input 
file (see Section 3.2.5, block 10). If the fan is rotating counter-clockwise viewed from the rear, then the 
rotational speed must be positive. Likewise, if the fan is rotating clockwise viewed from the rear, then the 
rotational speed must be negative. The program will check for this consistency, and will attempt to correct 
any errors. However, this check is not foolproof. The user should closely examine the velocity triangles 
that are displayed on the screen when Setup BFaNS is executed. 

The contents of the viscous CFD or measured data files can be viewed with a MATLAB (The 
MathWorks, Inc.) utility provided with Setup BFaNS. This utility consists of three files (reduce.m, 
stress.m, and only_plot.m) located in the directory that contains the Setup BFaNS source code. At the 
MATLAB prompt, type “reduce” and respond to the questions. 


TABLE 4.— FORMAT OF FLOW-FIELD FILES 
[Viscous CFD prediction or measured data] 


Block 

Variable Name 

Type 

Format 

Description 

1 

iwaketyp 

integer 

free 

not used 

2 

omega 

real 

free 

rotor rotational speed (rpm) 


xprtip 

real 

free 

axial distance from the rotor trailing edge 
tip to the velocity measurement plane tip 
location (inches) 

rhub 

real 

free 

inner duct radius at the velocity 
measurement plane (inches) 

rtip 

real 

free 

outer duct radius at the velocity 
measurement plane (inches) 

npassage 

integer 

free 

not used 


phipr 

real 

free 

not used 

3 

numr 

integer 

free 

number of radial measurement points for 
velocities 

numt 

integer 

free 

number of tangential measurement points for 
velocities 

4 

z (i) 

i=l , numr 

real 

free 

axial coordinates for velocities (inches) 
measured from blade tip trailing edge 

5 

x(i, j) 
i=l , numr 
j =1 , numt 

real 

free 

x-coordinates for velocities (inches) 

6 

y (i. j) 

i=l , numr 
j =1 , numt 

real 

free 

y-coordinates for velocities (inches) 

7 

icz 

integer 

free 

=1 if cz is in file, =0 if not 

icu 

integer 

free 

=1 if cu is in file, =0 if not 

icr 

integer 

free 

=1 if cr is in file, =0 if not 

itke 

integer 

free 

=1 if tke is in file, =0 if not 

ieps 

integer 

free 

=1 if eps is in file, =0 if not 

icpzcpz 

integer 

free 

=1 if cpzcpz is in file, =0 if not 

icpucpu 

integer 

free 

=1 if cpucpu is in file, =0 if not 

icprcpr 

integer 

free 

=1 if cprcpr is in file, =0 if not 

icpzcpu 

integer 

free 

=1 if cpzcpu is in file, =0 if not 

icpzcpr 

integer 

free 

=1 if cpzcpr is in file, =0 if not 

icpucpr 

integer 

free 

=1 if cpucpr is in file, =0 if not 

iscale 

integer 

free 

=1 if scale is in file, =0 if not 

8 

cz (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

axial component of velocity (ft/sec) 

9 

cu (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

circumferential component of velocity 
(ft/sec) 

10 

cr (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

radial component of velocity (ft/sec) 
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TABLE 4.— FORMAT OF FLOW-FIELD FILES 
[Viscous CFD prediction or measured data] 


Block 

Variable Name 

Type 

Format 

Description 

11 

tke (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

turbulent kinetic energy (ft**2/sec**2 ) 

12 

eps (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

turbulent dissipation rate (ft**2/sec**3) 

13 

cpzcpz (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

z-z reynolds stress (ft**2/sec**2) 

14 

cpucpu (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

0-0 reynolds stress (f t**2 /sec**2 ) 

15 

cprcpr (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

r-r reynolds stress (ft**2/sec**2 ) 

16 

cpzcpu (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

z-0 reynolds stress (ft**2/sec**2) 

17 

cpzcpr (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

z-r reynolds stress (ft**2/sec**2 ) 

18 

cpucpr (i, j ) 
i=l , numr 
j =1 , numt 

real 

free 

0-r reynolds stress (ft**2/sec**2 ) 

19 

scale (i) 
i=l , numr 

real 

free 

integral length scale (in) 

20 

ztip2 

real 

free 

axial distance from the rotor trailing edge 
tip to the PT2,TT2 measurement plane 
(inches) 

rhub2 

real 

free 

inner duct radius at the PT2,TT2 
measurement plane (inches) 

rtip2 

real 

free 

outer duct radius at the PT2,TT2 
measurement plane (inches) 

21 

numr2 

integer 

free 

number of radial measurement points for 
PT2 , TT2 

22 

r2 (i) 
i=l , numr2 

real 

free 

radial coordinate of radial measurement 
points for PT2 , TT2 (inches) 

23 

z2 (i) 
i=l , numr2 

real 

free 

axial coordinate relative to the rotor tip 
t.e. of radial measurement points for PT2, 
TT2 (inches) 

24 

pt2 (i) 
i=l , numr2 

real 

free 

circumferentially averaged total pressure 
(psfa) 

25 

tt2 (i) 
i=l , numr2 

real 

free 

circumferentially averaged total 
temperature (deg. R) 


TABLE 5.— FORMAT OF FLOW-FIELD FILES 


[Inviscid streamline prediction] 


Block 

Variable Name 

Type 

Format 

Description 

1 

f iledescr 

character 

a50 

file description 

2 

numr 

integer 

free 

number of radial measurement points 

3 

parmname 

character 

a50 

parameter name for block 4 

4 

diam (i) , 
i=l , numr 

real 

free 

local diameter (in) 

5 

parmname 

character 

a50 

parameter name for block 6 

6 

z (i) , 
i=l , numr 

real 

free 

axial coordinate (in) measured from flow 
path reference 

7 

parmname 

character 

a50 

parameter name for block 8 

8 

u(i) , 
i=l , numr 

real 

free 

local blade speed (in) 

9 

parmname 

character 

a50 

parameter name for block 10 

10 

cz (i) , 
i=l , numr 

real 

free 

axial component of velocity (ft/sec) 

11 

parmname 

character 

a50 

parameter name for block 12 
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TABLE 5.— FORMAT OF FLOW-FIELD FILES 


[Inviscid streamline prediction] 


Block 

Variable Name 

Type 

Format 

Description 

12 

Ct (i) , 
i=l , numr 

real 

free 

tangential component of velocity (ft/sec) 

13 

parmname 

character 

a50 

parameter name for block 14 

14 

cr (i) , 
i=l , numr 

real 

free 

radial component of velocity 

15 

parmname 

character 

a50 

parameter name for block 16 

16 

loss (i) , 
i=l , numr 

real 

free 

loss coefficient 

17 

parmname 

character 

a50 

parameter name for block 18 

18 

gapchord (i) , 
i=l , numr 

real 

free 

local gap/chord ratio 

19 

parmname 

character 

a50 

parameter name for block 20 

20 

pt2 (i) , 
i=l , numr 

real 

free 

total pressure (psfa) 

21 

parmname 

character 

a50 

parameter name for block 22 

22 

tt2 (i) , 
i=l , numr 

real 

free 

total temperature (deg. r) 


3.2.5 Create SetupBFaNS Input File 

The user must create a Setup BFaNS input file in the working sub-directory (i.e., the sub-directory 
where Setup BFaNS will be run). Table 6 explains the format of the Setup BFaNS input file, and 0 
shows a sample file. 

When zero is entered in blocks 21, 22, 23 or 24 of the input file, the corresponding boundary-layer 
thickness will be calculated automatically. For every new case, the author recommends that the user enter 
zeroes in each of these blocks. After running Setup BFaNS, the user should then confirm that the 
computed boundary layers are reasonable (see Section 3.2.8). If not, the boundary layers must be entered 
manually, and Setup BFaNS must be run again. See Section 4.10.1 for details on how Setup BFaNS 
estimates the boundary-layer thicknesses. 


TABLE 6.— FORMAT OF SETUP BFANS INPUT FILE 


Block 

Variable 

name 

Type 

Fonnat 

Description 

1 

filenam 

character 

a30 

Specifies the name of the file that contains the 
flow path coordinates. The length of the filename 
must be less than or equal to 30 characters. This 
file must be located in the directory specified by 
Card 5002 of the Setup BFaNS author file (note: the 
default geometry directory is bfans_geometry) . 

2 

zblade 

real 

free 

Specifies the axial coordinate (in inches) of the 
blade trailing edge tip. This coordinate must be 
measured from the same reference point as the axial 
coordinates of the flow path. 

3 

zvane 

real 

free 

Specifies the axial coordinate (in inches) of the 
vane trailing edge tip. This coordinate must be 
measured from the same reference point as the axial 
coordinates of the flow path. 

4 

zhigh 

real 

free 

Specifies the axial coordinate (in inches) of the 
hilite. This coordinate must be measured from the 
same reference point as the axial coordinates of 
the flow path. 

5 

znose 

real 

free 

Specifies the axial coordinate (in inches) of the 
nose cone leading edge. This coordinate must be 
measured from the same reference point as the axial 
coordinates of the flow path. 

6 

zcore 

real 

free 

Specifies the axial coordinate (in inches) of the 
splitter leading edge. This coordinate must be 
measured from the same reference point as the axial 
coordinates of the flow path. 
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TABLE 6.— FORMAT OF SETUP BFANS INPUT FILE 


Block 

Variable 

name 

Type 

Format 

Description 

7 

zexit 

real 

free 

Specifies the axial coordinate (in inches) of the 
duct exit. This coordinate must be measured from 
the same reference point as the axial coordinates 
of the flow path. 

8 

filenam 

character 

a30 

Specifies the name of the file that contains the 
blade geometry. The length of the filename must be 
less than or equal to 30 characters. This file 
must be located in the directory specified by Card 
5002 of the Setup BFaNS author file (note: the 
default geometry directory is bfansgeometry) . 

9 

bnum 

integer 

free 

Specifies the number of fan blades. 

10 

dir 

integer 

free 

Specifies the fan's direction of rotation. Enter 
+1 if the fan is rotating counter-clockwise viewed 
from the rear, or -1 if the fan is rotating 
clockwise viewed from the rear. This notation is 
consistent with the cylindrical coordinate system 
shown earlier, where 0 increases in the counter- 
clockwise direction viewed from the rear of the 
engine . 

11 

filenam 

character 

a30 

Specifies the name of the file that contains the 
vane geometry. The length of the filename must be 
less than or equal to 30 characters. This file 
must be located in the directory specified by Card 
5002 of the Setup BFaNS author file (note: the 
default geometry directory is bfans_geometry) . 

12 

vnum 

integer 

free 

Specifies the number of vanes. 

13 

filenam 

character 

a30 

Specifies the name of the file that contains the 
flow-field description at Station #1 upstream from 
the fan. The length of the filename must be less 
than or equal to 30 characters. 

14 

filetypl 

integer 

free 

Specifies whether the file specified in Block 13 
contains information from an inviscid streamline 
deck, a viscous cfd deck, or measured data. Enter 
1 for streamline, 2 for cfd, 3 for data. If 1, the 
file specified by Block 13 must be located in the 
directory specified by Card 5004 of the Setup BFaNS 
author file. If 2, the file specified by Block 13 
must be located in the directory specified by Card 
5006 of the Setup BFaNS author file. If 3, the 
file specified by Block 13 must be located in the 
directory specified by Card 5008 of the Setup BFaNS 
author file. 

15 

filenam 

character 

a30 

Specifies the name of the file that contains the 
flow-field description at Station #2 downstream 
from the fan. The length of the filename must be 
less than or equal to 30 characters. 

16 

filetyp2 

integer 

free 

Specifies whether the file specified in Block 15 
contains information from an inviscid streamline 
deck, a viscous cfd deck, or measured data. Enter 
1 for streamline, 2 for cfd, 3 for data. If 1, the 
file specified by Block 15 must be located in the 
directory specified by Card 5004 of the Setup BFaNS 
author file. If 2, the file specified by Block 15 
must be located in the directory specified by Card 
5006 of the Setup BFaNS author file. If 3, the 
file specified by Block 15 must be located in the 
directory specified by Card 5008 of the Setup BFaNS 
author file. 

17 

filenam 

character 

a30 

Specifies the name of the file that contains the 
flow-field description at Station #3 upstream from 
the vane. The length of the filename must be less 
than or equal to 30 characters. 
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TABLE 6.— FORMAT OF SETUP BFANS INPUT FILE 


Block 

Variable 

name 

Type 

Fonnat 

Description 

18 

filetyp3 

integer 

free 

Specifies whether the file specified in Block 17 
contains information from an inviscid streamline 
deck, a viscous cfd deck, or measured data. Enter 
1 for streamline, 2 for cfd, 3 for data. If 1, the 
file specified by Block 17 must be located in the 
directory specified by Card 5004 of the Setup BFaNS 
author file. If 2, the file specified by Block 17 
must be located in the directory specified by Card 
5006 of the Setup BFaNS author file. If 3, the 
file specified by Block 17 must be located in the 
directory specified by Card 5008 of the Setup BFaNS 
author file. 

19 

filenam 

character 

a30 

Specifies the name of the file that contains the 
flow-field description at Station #4 downstream 
from the vane. The length of the filename must be 
less than or equal to 30 characters. 

20 

filetyp4 

integer 

free 

Specifies whether the file specified in Block 19 
contains information from an inviscid streamline 
deck, a viscous cfd deck, or measured data. Enter 
1 for streamline, 2 for cfd, 3 for data. If 1, the 
file specified by Block 19 must be located in the 
directory specified by Card 5004 of the Setup BFaNS 
author file. If 2, the file specified by Block 19 
must be located in the directory specified by Card 
5006 of the Setup BFaNS author file. If 3, the 
file specified by Block 19 must be located in the 
directory specified by Card 5008 of the Setup BFaNS 
author file. 

21 

del 

real 

free 

Specifies the boundary layer thickness (in inches) 
on the shroud surface at Station #1. If zero is 
entered, the boundary layer thickness will be 
calculated automatically 

22 

del 

real 

free 

Specifies the boundary layer thickness (in inches) 
on the hub surface at Station #1. If zero is 
entered, the boundary layer thickness will be 
calculated automatically 

23 

del 

real 

free 

Specifies the boundary layer thickness (in inches) 
on the shroud surface at Station #3. If zero is 
entered, the boundary layer thickness will be 
calculated automatically 

24 

del 

real 

free 

Specifies the boundary layer thickness (in inches) 
on the upper splitter surface at Station #3. If 
zero is entered, the boundary layer thickness will 
be calculated automatically 

25 

titlel 

character 

a30 

Specifies the first line of the run title (must be 
less than or equal to 30 characters) 

26 

title2 

character 

a30 

Specifies the second line of the run title (must be 
less than or equal to 30 characters) 


3.2.6 Create Setup BFaNS Author File (Optional) 

The author file allows the user to change default values that are hardcoded in Setup BFaNS. 

Setup BFaNS will search the working directory for a file named “setupbfans. author”. If this file exists, 
Setup BFaNS will use it to over-ride the internal defaults for various options. Appendix J shows a sample 
Setup BFaNS author file. The MATLAB program “sbf author.m” creates an author-file template. This 
program is located in the Setup BFaNS source code directory. 

The author file has distinct blocks of input referred to as cards, with each card beginning on a new 
line. A negative card number indicates that the card is turned off. A positive card number indicates that a 
card is turned on. This feature allows the user to activate/deactivate various options in Setup BFaNS. 
Each card has five input parameters (two integer, two real and one alphanumeric), which are provided at 
the end of the line. When a card is turned on, its input parameters are accessible to every subroutine 
within Setup BFaNS. 
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TABLE 7.— FORMAT OF SETUP BFANS AUTHOR FILE 


Card 

Description 

5000 

If turned on, this card specifies the file name of the Setup BFaNS input file, and 
the user will not be prompted to enter the file name. 

5002 

If turned on, this card specifies the name of the sub-directory that contains the 
geometry files (flow path, blade, and vane) . The default sub-directory is 

bf ans_geometry . 

5004 

If turned on, this card specifies the name of the sub-directory that contains the 
flow-field files generated with an inviscid streamline deck. The default sub- 
directory is bfans sline . 

5006 

If turned on, this card specifies the name of the sub-directory that contains the 
flow-field files generated with a viscous cfd deck. The default sub-directory is 

bfans cfd. 

5008 

If turned on, this card specifies the name of the sub-directory that contains the 
flow-field files generated with measured data. The default sub-directory is 

bf ans_data . 

5010 

If turned on, this card specifies the path to the Setup BFaNS source code. The 
default path is ~/Codes/fans_nasa/setup_bfans . 

5012 

through 

5048 

Reserved for future use 

5050 

If turned on, a debug file (debug.dat) will be created in the working directory. 

5052 

If turned on, the turbulent length scale at Station #1 will be determined in the 
following manner: 

1. if imoda(5052) = 1, the turbulent length scale will be computed from the k- 
epsilon turbulence model (see cards 5056, 5058, and 5060) . 

2. if imoda(5052) = 2, the turbulent length scale will be determined from 
empirical correlations (see cards 5062 and 5064) that relate scale to 
boundary layer thickness. 

The default is to use the turbulent length scales specified in the flow-field input 
file . 

5054 

If turned on, the turbulent length scale at Station #3 will be determined in the 
following manner: 

1. if imoda(5054) = 1, the turbulent length scale will be computed from the k- 
epsilon turbulence model (see cards 5056, 5058, and 5060) . 

2. if imoda(5054) = 2, the turbulent length scale will be determined from 
empirical correlation (see cards 5066, 5068, and 5070) that relate scale to 
boundary layer thickness and wake thickness. 

The default is to use the turbulent length scales specified in the flow-field input 
file . 

5056 

If turned on, this card specifies the coefficient that relates the turbulence 
length scale to the k-epsilon model. The coefficient is specified by the third 
parameter (real) . The default coefficient is one. This card is used only when Card 
5052 is on, and imod(5052) = 1. 

5058 

If turned on, this card will multiply the turbulent kinetic energy (K) by the ratio 
specified by the third parameter (real) . This card only applies to viscous CFD 
input. This card is used only when Card 5052 is on, and imod(5052) = 1. 

5060 

If turned on, this card will multiply the turbulent dissipation (epsilon) by the 
ratio specified by the third parameter (real) . This card only applies to viscous 
CFD input. This card is used only when Card 5052 is on, and imod(5052) = 1. 

5062 

If turned on, this card specifies the ratio between the turbulent length scale and 
the hub boundary layer thickness at Station #1. The default ratio is 0.62. This 
card is used only when Card 5052 is on, and imod(5052) = 2. 

5064 

If turned on, this card specifies the ratio between the turbulent length scale and 
the shroud boundary layer thickness at Station #1. The default ratio is 0.62. 

This card is used only when Card 5052 is on, and imod(5052) = 2. 

5066 

If turned on, this card specifies the ratio between the turbulent length scale and 
the upper splitter boundary layer thickness at Station #3. The default ratio is 
0.62. This card is used only when Card 5054 is off, and imod(5054) = 2. 

5068 

If turned on, this card specifies the ratio between the turbulent length scale and 
the shroud boundary layer thickness at Station #3. The default ratio is 0.62. 

This card is used only when Card 5054 is off, and imod(5054) = 2. 

5070 

If turned on, this card specifies the ratio between the turbulent length scale and 
the wake thickness at Station #3. The default ratio is 0.68. This card is used 
only when Card 5054 is off, and imod(5054) = 2. 

5072 

through 

5078 

Always turn these off. These cards can be used to curve fit the wake centerline. 

5080 

If turned on, streamtube contraction effects will be included in the calculation of 
cascade parameters 
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TABLE 7.— FORMAT OF SETUP BFANS AUTHOR FILE 


Card 

Description 

5082 

If turned on, upwash will be computed from the turbulent kinetic energy assuming 
isotropic turbulence 

5084 

If turned on. Setup BFaNS will not check to see if cfd or experimental 
covers one full pitch 

9999 

Indicates end of author file. 


3.2.7 Execute SetupBFaNS 

Table 8 shows the procedure for executing Setup BFaNS. 


TABLE 8.— EXECUTING SETUP BFANS 


Step 1 

Follow the instruction in Sections 3.2.1 through 3.2.6. 

Step 2 

Go to the working directory. 

Step 3 

Execute Setup BFaNS using the alias sbf. 

Step 4 

When prompted, enter the name of the Setup BFaNS author file. 

Step 5 

When prompted, enter the name of the Setup BFaNS input file (unless Card 5000 in 
the Setup BFaNS author file is turned on, in which case the file name is specified 
in the author file) . 

Step 6 

A figure will be displayed on the screen showing the blade velocity triangles. 
Look at the figure and confirm that the velocity triangles are correct relative to 
the blade orientation. If satisfied, place your cursor inside the figure window, 
and type the letter x on your keyboard. Typing x will continue program execution. 

Step 7 

A figure will be displayed on the screen showing the vane velocity triangles. Look 
at the figure and confirm that the velocity triangles are correct relative to the 
vane orientation. If satisfied, place your cursor inside the figure window, and 
type the letter x on your keyboard. Typing x will continue program execution. 

Step 8 

A figure will be displayed on the screen showing the fan-duct configuration. Look 
at the figure and confirm that the configuration is correct. If satisfied, place 
your cursor inside the figure window, and type the letter x on your keyboard. 
Typing x will continue program execution. 

Step 9 

A menu will be displayed to allow you to change the geometry. If you wish to 
change something, type in the number corresponding to your selection and press 
return. Then provide the requested input. The old and new geometry will now be 
shown on the screen. Place your cursor inside the figure window, and type the 
letter x on your keyboard. Typing x will continue program execution, allowing you 
to make more changes to the geometry. 

Step 10 

After all geometric changes have been completed, type 99 to finish running 
Setup BFANS. 

Step 11 

Confirm that you have a file named bfans. input in the working directory, along with 
a sub-directory . 


3.2.8 Examine Boundary-Layer Characteristics 

For every new case, the author recommends that the user enter zeroes in blocks 21, 22, 23 and 24 of 
the input file. By entering zeroes, the endwall boundary-layer thicknesses will be calculated 
automatically. However, the automated calculation procedure is not very robust, so the user should 
confirm that the computed boundary layers are reasonable. If not, the boundary layers must be entered 
manually, and Setup BFaNS must be run again. 

Table 9 shows the procedure for checking the endwall boundary layers. This procedure assumes that 
the user has access to MATLAB. 


TABLE 9.— EXAMINING ENDWALL BOUNDARY LAYERS 


Step 1 

Go to the MATLAB sub-directory inside the working sub-directory. This sub-directory 
should contain several m-files, including "plot bl.m" , "stationl .m" , and 
"station3 ,m" . 

Step 2 

Launch MATLAB 

Step 3 

At the MATLAB prompt, type plot bl.m 

Step 4 

The boundary-layer profiles at Station #1 and Station #3 will be displayed. If the 
boundary layers look incorrect, visually estimate the correct boundary layer 
thickness . 

Step 5 

Enter the correct boundary-layer thicknesses in the Setup BFaNS input file (Blocks 
21 - 24), and re-run Setup BFaNS . 
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3.2.9 Examine Wake Characteristics 


Table 10 shows the procedure for checking the fan wake characteristics. This procedure assumes that 
the user has access to MATLAB. This step is important only when the turbulent length scales at Station 
no. 3 are being computed from the fan wake thickness (as determined by Card 5054 in the author file). 


TABLE 10.— EXAMINING WAKE CHARACTERISTICS 


Step 1 

Go to the MATLAB sub-directory inside the working directory. This sub-directory 
should contain several m-files, including "plot wake.m". 

Step 2 

Launch MATLAB 

Step 3 

At the MATLAB prompt, type plot wake.m. 

Step 4 

Examine the wake profiles at each radius, focusing on the region away from the 
endwalls . 

Step 5 

If the wakes look incorrect, contact the author. Currently, there is no way to 
manually enter the wake thicknesses. 


3.2.10 Examine Output File 

Table 1 1 provides suggestions for examining the SetupBFaNS output file. These suggestions are 
based on some of the problems that the author encountered when first running Setup BFaNS. These 
problems usually indicated an error in one of the flow-field input files (e.g., in the way that the flow-field 
predictions were performed). 


TABLE 1 1.— EXAMINING THE OUTPUT FILE 


Step 1 

Confirm that there are no overflows (i.e *****) in the output file 

Step 2 

Confirm that there are no reverse flow regions (Cz < 0) 

Step 3 

Confirm that the swirl is zero at the rotor leading edge (unless an inlet guide 
vane is present) 

Step 4 

Confirm that the swirl between the rotor and stator has the same sign as the 
rotational speed of the rotor 

Step 5 

Confirm that there is little or no swirl at the vane trailing edge 

Step 6 

Confirm that the swirl at the hub of the fan leading edge is equal to the hub 
rotation speed 

Step 7 

Confirm that the turbulent scales in the endwall regions are less than the 
boundary-layer thickness; confirm that the turbulent length scales in the wake 
region are less than the wake thickness 


4.0 Program Documentation 

This section provides technical documentation for the Setup BFaNS computer program, which is 
written in FORTRAN 77. Appendix L shows the subroutine hierarchy. 

4.1 Units 

Program inputs and outputs are in English units. However, internal calculations are done in SI units. 

4.2 Geometry Input 

Geometric information is passed into the main program through subroutine “get geom”. This routine 
first reads the flow-path geometry, then the blade geometry, and finally the vane geometry. 

4.2.1 Flow-Path Geometry 

The flow-path geometry is read by the subroutine “get duct geom”. This routine reads the (z,r) 
coordinates that define the hub surface, the shroud surface, the lower splitter surface and the upper splitter 
surface. This routine also reads the axial location of the blade tip trailing edge (ZBLADE), the vane tip 
trailing edge (ZVANE), the hilite (ZH1GH), the nose-cone leading edge (ZNOSE), the splitter leading 
edge (ZCORE), and the fan-duct trailing edge (ZEX1T). See Figure 5 for a description of the flow-path 
terminology. 
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4.2.2 Blade and Vane Geometry 

The blade and vane geometries are read by the routine “get row geom”, which is called separately 
for the blade and the vane. This routine reads the (x*,y*,z*) coordinates and metal angle of the leading 
and trailing edges of the airfoil. See Figure 6 for a description of the airfoil coordinate system. 

After reading the airfoil geometry, “get row geom” applies the following transformations: 

• Rotates the airfoil to ensure that the innermost point on the leading edge is located on the +x axis, 


9/£i = tan 


-l 


r * \ 
yiD 


\ X ID 

* r\* * • rv* 

x = x cos vjD+y sint ) jp 
y — -x sin vjD+y cos 0 jp 


Converts from Cartesian coordinates to cylindrical coordinates, i.e., 

r - tJx 2 + y 2 

tan 0 = — 
x 

• Ensures that the airfoil is oriented properly relative to the direction of rotation (this process is 
referred to as “flipping” the airfoil). “Flipping” is required because SetupBFaNS internally uses 
metal angles that are measured from the +0 direction, whereas the airfoils are defined relative to 
their suction-side tangential. 

• If the fan is rotating counter-clockwise (as determined by Block 10 in the Setup BFaNS input 
file), Setup BFaNS will take the supplement of the blade metal angles. This process is needed 
because the suction-side tangential of the blade is in the -0 direction for counter-clockwise 
rotation. In addition, Setup BFaNS will multiply the 0-coordinates of the blade by — 1. No 
changes are made to the vane geometry because its suction-side tangential is already in the +0 
direction. 

• If the fan is rotating clockwise, Setup BFaNS will take the supplement of the vane metal angles. 
This process is needed because the suction-side tangential of the vane is in the -0 direction for 
clockwise rotation. In addition, Setup BFaNS will multiply the 0-coordinates of the vane by — 1. 
No changes are made to the blade geometry because its suction-side tangential is already in the 
+0 direction. 

• Translates the blade and vane to ensure that their tip trailing edges are located at ZBLADE and 
ZVANE, respectively, 


z - z* + ZBLADE or z = z* + ZVANE 


4.3 Flow-Field Input 

Flow-field information is passed into the main program through subroutine “get flow data”. This 
routine reads the velocity, total-pressure and total-temperature fields at four stations (see Figure 5): 

• Station 1 : upstream from the blade leading edge 

• Station 2: downstream from the blade trailing edge (but upstream from the splitter leading edge) 

• Station 3: upstream from the vane leading edge (but downstream from the splitter leading edge) 

• Station 4: downstream from the vane trailing edge. 
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It also circumferentially averages the flow field, computes the density and speed-of-sound profiles, 
and computes the streamfunction. Finally, “get flow data” ensures that the fan speed is in the direction 
specified in the Setup BFaNS input file. 

4.3.1 Velocity Field 

The coordinates of the velocity grid are read by the subroutine “get mesh data”. This routine reads 
the (x*,y*,z*) coordinates at which velocity data is available (the velocity grid must define an 
axisymmetric surface). When using experimental data or viscous CFD, the z*-coordinates must be 
measured from the blade tip trailing edge. When using results from an inviscid streamline deck, the z*- 
coordinates must be measured from the same reference that was used to define the flow-path geometry. 

After reading the data points, “get mesh data” applies the following transformations: 


• Rotates the grid to ensure that the first point is located on the +x axis. This transformation is 
permissible because all flow -field quantities will ultimately be averaged in the circumferential 
direction (i.e., the location of the wake relative to the blade or vane is not important for 
broadband noise). 


0/0 = tan 


-1 


/ * \ 

yip 

V X ID y 


* * • /-v* 

x = x cos ujd + y sin 0 id 
y — -x sm 0/o + y cos 0 /o 


• Converts from Cartesian coordinates to cylindrical coordinates, i.e., 


r = -Jx 2 +y 2 

tan 0 = — 

x 


• If necessary, translates the grid to ensure that the tip is located at the correct axial location, i.e., 

z = z* + ZBLADE 

This translation is not done when the flow- field is provided by an inviscid streamline prediction. 

Once the grid coordinates are known, the mean and turbulent velocity fields are read by the 
subroutine “get velocity data”. 

When the flow-field is based on experimental data or viscous CFD, and some components of the 
Reynolds-stress tensor are not available, “get velocity data” estimates them based on the turbulence 
information that is available (see Table 12). The available turbulence information must satisfy one of the 
following conditions: 


Condition 

Known turbulence quantities 

1 

Turbulent kinetic energy 

2 

One nonnal-stress component 

3 

Two nonnal-stress components 

4 

Two nonnal-stress components, one shear-stress component 

5 

Three normal-stress components 

6 

Three normal-stress components, one shear-stress component 

7 

Three normal-stress components, two shear-stress components 

8 

Three normal-stress components, three shear-stress components 
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These conditions are consistent with the turbulence data that might be available from hot-wire or laser 
Doppler velocimetry (LDV) experiments. 

If necessary, “get velocity data” extends the grid to the walls, where it sets the velocity to zero 
(relative to the wall), along with all turbulence quantities. For viscous-CFD results, “get velocity data” 
also checks for negative values of K (which are physically impossible) and sets them and the 
corresponding Reynolds-stress components equal to zero. This step eliminates numerical problems later 
in the computer program. 

For viscous-CFD results, the user has the option to multiply K and all components of the Reynolds- 
stress tensor by the constant specified in Card 5058 of the author file. The user also has the option to 
multiply the viscous dissipation (s) by the constant specified in Card 5060 of the author file. These 
options can be used to “adjust” the turbulence predictions to “improve” agreement with measured data. 

4.3.2 Circumferential Averaging 

Subroutine “circum average” circumferentially averages the velocity field in the following manner, 

. 9 max 

Qavg = a Tfj J 

t'max ^min n 

^min 

In the equation above, Q can represent any velocity component ( C,j, any Reynolds-stress component 
( cfij ), turbulent kinetic energy (AT), or viscous dissipation (s). The integration uses the trapezoidal rule. 


TABLE 12.— ESTIMATING THE REYNOLDS-STRESS TENSOR 


Condition 

K 

c z c z 

c e c e 

c r c r 

c z c e 

C-C r 

c B c r 

1 

known 

= 2 -k 
3 

= !* 
3 

3 

0 

0 

0 

2 

3 

= — C 7 C 7 
2 

known 

= C Z C z 

= C Z C z 

0 

0 

0 


3 

= — c e c 0 

= c e c e 

known 

= c e c e 

0 

0 

0 


3 

= —c r c, 
2 

= c r c r 

= CfCy 

known 

0 

0 

0 

3 

II 

+ 

known 

known 

= —( c z c z +cece) 

0 

0 

0 


3/ \ 

= — \C 7 C 7 + c r c r I 

4 ' ’ 

known 

= \(c.c z +c r c r ) 

known 

0 

0 

0 


^(ce^e +c r c r ) 

= —(cb c q +c r c r ) 

known 

known 

0 

0 

0 

4 

= 4tzC z +C e C 0 ) 

4 v 

known 

known 

= —( c z c z +C0C©) 

known 

= C z Cq 

= C z Cq 


3/ 1 

= — 1 c,c, + c r c r 1 
4 V " ’ 

known 

= j(czC z +c r c r ) 

known 

= C z C r 

known 

= c z c r 


= — (cece + <v<v) 

= ^(cece +c r c r ) 

known 

known 

= c e c r 

= C Q C r 

known 

5 

= —{c z C z + CqCq C r C r j 

known 

known 

known 

0 

0 

0 

6 

= j (c z C z + CqCq + C r C r ) 

known 

known 

known 

known 

= C z Cq 

= C z Cq 


= — (c z c z + CqCq + C r C r j 

known 

known 

known 

= c z c r 

known 

= c z c r 


= ~^[c z C z + CqCq + C r C r ) 

known 

known 

known 

= c e c r 

= CQC r 

known 

7 

= — (c z C z + CqCq -\-C r C r j 

known 

known 

known 

known 

known 

= \(c z c e +c ! c r ) 


= y (c z C z + CqCq + C r C r ) 

known 

known 

known 

known 

= —{c z Cq +C0C r ) 

known 


= ~( c z c z T CqCq + C r C r j 

known 

known 

known 

= j(c z c r +c e c r ) 

known 

known 

8 

= — (c z C z + CqCq + C r C r ) 

known 

known 

known 

known 

known 

known 
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4.3.3 Computing Density and Speed of Sound 

Density and speed of sound are computed in subroutine “getrhosound”. “Getrhosound” reads the 
total-pressure and total-temperature profiles at Stations 1 to 4, and computes density and speed of sound 
via: 


T = T t 

P = P T 

P = 
a = 


1 cl 


avg 


2 C p 

(T t ) 

-y 
| y — i 

[t) 


P 


RT 



Vr RT 


The specific heat (C p ), gas constant (R) and specific-heat ratio (y)are hardcoded in “get rho sound” to be 
1005 m 2 /s 2 -K, 287 m 2 /s 2 -K and 1.4, respectively (i.e., the values for air — see Ref. 8). 

In some instances, the total-pressure, total-temperature and velocity data are acquired at different 
axial and radial locations. To account for these instances, “get rho sound” interpolates the total-pressure 
and total-temperature data onto radii that correspond to the same area fractions that the velocity data were 
acquired, i.e., 


( 




interp ID 


r OD r ID 




J 


Pf,Tj Data 


( r 2 - r 2 ^ 
measured ID 


r, 


OD r ID 


'Velocity Data 


The values of r interp are computed, and the (Pr, 7V) data are interpolated onto those radii. Figure 8 
demonstrates this process, as it would apply to Station no. 3. The total pressure and total temperature are 
assumed to be constant along lines of constant area fraction (i.e., the lines connecting the x’s and o’s). 
The assumption is that lines of constant area fraction approximate the streamtubes. 

Station #3 



N ASA/CR— 20 1 0-2 1 6898/VOL 1 


20 



4.3.4 Computing Streamfunction 

Streamfunction is computed in subroutine “stream function” by integrating the mass flux profile 
from the outer wall to the inner wall via: 


S(r)= Jp Cayg-dA 

? od 

Provided there is no reverse flow, there will be a one-to-one correspondence between radius and 
streamfunction. 

Note that the value of the streamfunction at the ID corresponds to the mass flow rate through the 
measurement surface. 

4.3.5 Correcting Direction of Rotation 

“Getflowdata” determines whether the fan speed specified in the flow-field input file is consistent 
with the direction specified in the SetupBFaNS input file. If the fan speed is in the wrong direction, then 
“get flow data” makes the following corrections: 


Q = -Q 

Q = -c e 


c z CQ — C z Cq 
CQCf — —CQCf 


These corrections are also made to the corresponding circumferentially averaged quantities. 

4.4 Checking Consistency of Fan Speeds 

After reading the flow-field data at each station, Setup BFaNS checks the consistency between the 
fan speeds that were read. If the fan speed at any station is different by 1 to 2 percent of the average fan 
speed, the program issues a warning and continues. If the fan speed at any station is different by more 
than 2 percent of the average fan speed, the program is terminated. 

4.5 Modifications to Flow Path 

To facilitate interpolation of the flow field, the subroutine “add_points” makes the following changes 
to the flow path: 

• Shroud Surface . — Adds the outside-diameter coordinates of the blade leading edge and trailing 
edge; adds the outside-diameter coordinates of the vane leading edge and trailing edge; adds the 
outside-diameter coordinates of the velocity grid at Stations 1 to 4 

• Hub Surface . — Adds the inside-diameter coordinates of the blade leading edge and trailing edge; 
adds the inside-diameter coordinates of the velocity grid at Stations 1 and 2 

• Upper Splitter Surface . — Adds the inside-diameter coordinates of the vane leading edge and 
trailing edge; adds the inside-diameter coordinates of the velocity grid at Stations 3 and 4 

4.6 Streamline Estimation 

The blade geometry, vane geometry and flow field are mapped onto lines of constant streamfunction 
(i.e., streamtubes). This mapping is done in four steps: 
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• Subroutine “normalize” normalizes the streamfunction by the mass flow rate passing through the 
vane. In this manner, the normalized streamfunction varies between 0 and 1 at Stations 3 and 4, 
and between 0 and 1+BPR at Stations 1 and 2 (where BPR = bypass ratio). 

If the flow rate through Stations 1 and 2 are different by more than 1 percent, this routine warns 
the user. Likewise, this routine warns the user if the flow rate through Stations 3 and 4 are 
different by more than 1 percent. It also warns the user if there is more flow passing through the 
stator than the rotor (i.e., when BPR < 0). These warnings indicate an error in the input files, 
though the program will continue to run. 

• Subroutine “stream line” maps Station 2 flow-field information onto the streamtubes 
corresponding to Station 1 . The mapping is done by finding the radii at Station 2 that correspond 
to the normalized streamfunction values at Station 1 . A similar mapping is done between Stations 
4 and 3. This routine will fail if a one-to-one correspondence between radius and streamfunction 
does not exist (e.g., when there is reverse flow). 

• Subroutine “intersect2” estimates where the streamtubes intersect the blade and vane. It does this 
by determining the radii at the edge that correspond to the area fractions at the measurement 
station: 


( 


V 


.2 _ .2 
r interp r ID 


r OD r I~D 




J 


Edge 


( r 2 -r 2 ^ \ 
measured ID 


r OD r ID 


J Measurement 
Station 


The assumption is that lines of constant area fraction approximate the streamtubes. The axial 
coordinate, tangential coordinate and metal angle of the airfoil are interpolated onto these 
interpolated radii. 

• Subroutine “estimate2” estimates the flow-field information at the leading edge and trailing edge 
of the fan and FEGV. The only quantities affected are the axial, radial and tangential velocities. 
The axial and radial velocities are adjusted based on the change in area between the measurement 
station and the airfoil edge. The tangential velocity is adjusted based on the change in radius, 


Cr,edge ~ C r> station 


C%,edge Gt, station 


A edge J 
\ 

f station 
r edge y 


4.7 Confirm Blade and Vane Orientation 

Subroutine “triangle” displays the airfoil geometries and the corresponding inlet and exit velocity 
triangles on the screen. The user must examine the figures to confirm that the airfoils are oriented 
properly relative to the velocity triangles. If acceptable, the user must place the cursor in the figure 
window and hit “x” to continue. 

4.8 Airfoil Geometry Calculations 

Chord, gap, solidity, stagger angle, geometric sweep and geometric lean are computed in subroutine 
“geom”. The chord, solidity and stagger angle are computed on the “unwrapped” cascade along the 
(m,r0) surfaces. 
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4.8.1 Airfoil Chord 

Chord is defined by: 


b = V (x,e ~ X le ) 2 + [y,e ~ Vie f + [ Z ,e ~ Z le f 

Meridional chord is defined by: 

K = Vfc-'/J + ( Z te~ Z leY 

4.8.2 Airfoil Gap 

Gap is defined by: 

2nri 

x = — where N = number of airfoils 

N 


4.8.3 Solidity 

Cascade solidity is defined by: 


b 



Gap-to-chord ratio is equal to — . 

CT 

4.8.4 Stagger Angle 

Stagger angle is defined by: 


c, = tan 


-l 


Ste^te-Vety, 


e J 


In Setup BFaNS, the stagger angle is measured from the 0-direction. This definition is not consistent 
with Glegg or Hanson, who measure stagger angle from the axial direction. Therefore, when applying 
Hanson or Glegg’ s theories, we must replace their stagger angle with the complement of q. 

4.8.5 Geometric Sweep 

To compute geometric sweep at the fth point on the leading edge, Setup BFaNS determines the 
parabola (z as a function of r) that satisfies the following conditions: 

z = C2>' 2 +Cf' + Co 
z = z le,i~ 1 at r = r le,i~ 1 
z = z te>i at r = r tei 

z = z le,i+\ at r - ri e i+ i 

The coefficients are determined in subroutine “parabola”. From this parabola, the geometric sweep at 
the fth point is determined via, 
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dz 

tanv) >s= — 
dr 


The sweep at the first point is assumed to be equal to that at the second point. Likewise, the sweep at 
the last point is assumed to be equal to that at the second-to-last point. 

4.8.6 Geometric Lean 

To compute geometric lean at the f'th point on the leading edge, Setup BFaNS determines the parabola (9 
as a function of In r) that satisfies the following conditions: 


0 = i; 2 (lnr) 2 

+ Cl (in r)+ Co 

7 

c£T 

II 

cr> 

at 

lnr = in r lei _\ 

-sF 

c£r 

ii 

cr> 

at 

In r = In r leJ 

+ 

-$F 

CD 

II 

CD 

at 

In r = In r leJ+] 


The coefficients are determined in subroutine “parabola”. From this parabola, the geometric lean at the 
z'th point is determined via, 


tan\|/£ = 


1 dr 
r dO 


d (In r) 
dO 


The lean at the first point is assumed to be equal to that at the second point. Likewise, the lean at the 
last point is assumed to be equal to that at the second-to-last point. 


4.9 Flow-Field Calculations 

Leading-edge Mach number, meridional flow angle and incidence are computed in subroutine “flow”. 
Cascade parameters such as the lift, drag and loss coefficients are computed in subroutine “cascade”. 
Equivalent “isolated airfoil” angle-of-attack is computed in subroutine “naca0012”. 


4.9.1 Leading-Edge Mach Number 

Leading-edge Mach number is defined by: 




W, 


avg 


le 


for the rotor 


M k 


C, 


avg 


le 


for the stator 


4.9.2 Meridional Flow Angle 

The meridional flow angle is defined as the angle between a streamtube and the radial direction. Since 
the points that define the airfoil have been mapped onto streamtubes, then the meridional flow angle can 
be determined from: 


f . 


<|) = tan 1 


z le z /e 

l r te -n e J 
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4.9.3 Incidence 


Incidence is defined as the difference between the leading-edge flow angle and the leading-edge 
metal angle, both measured from the meridional direction. The leading edge flow angle is, 


P/e = tan 1 


\J^m,avg J 


le 


For a blade, the incidence is, 


For a vane, the incidence is, 


i = P te 



P/e 


A 


J 




P/e 


In both cases, the n/2 term is needed because the metal-angle (P*) is referenced to the 0-direction rather 
than the meridional direction. 


4.9.4 Cascade Parameters 

The cascade loss coefficient is determined from a correlation based on diffusion factor (Ref. 9 and 
10). The diffusion factor for a two-dimensional cascade is defined as, 


W te | IFg /e — Wft jg 

W le 2oW le 


(see Ref. 10, Eq. (54)) 


When there is a change in radius from the leading edge the trailing edge, the more general form for the 
diffusion factor is, 


D _ 1 K , r ieW eJ e- r t eWe,,e 

' a (r„+r,,)F„ 

Reference 1 0 presents a plot of wake momentum thickness (normalized by chord) versus diffusion factor 
for several sets of cascade data. A reasonably good fit to the data is provided by, 

= 0.006 + 0.0002(e 7 5jD / - 1) (see Ref. 10, Fig. 148) 

For a two-dimensional cascade with incompressible flow, the cascade loss, drag and lift coefficients can 
be determined from, 


co = 


2 

f s 2 \ 

f-ll 

( COS P/e N 

COS P/e 

U J 

UJ 

(cOSp te J 


^ _ x cos 3 p„, 

C D = ™T — TiT - 

b COS 7 P/ e 


(see Ref. 11, Eq. (3.32)) 
(see Ref. 11, Eq. (3.33)) 
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(see Ref. 11, Eq. (3.18)) 


C L =2^(tan(3 /e -tan(3 te )cos|3 w -Q,tan|3 m 
b 

where the mean cascade angle is defined by, 

= tan -1 { tan + tan 1 (see Ref. 1 1, Eq. (3.6)) 

V 2 

To avoid a divide-by-zero error in D t , these calculations are not done at the walls. Rather, the loss, lift 
and drag coefficients are assumed to be constant over the first two and last two points on the airfoil. 

These incompressible cascade relationships assume that the axial velocity at the trailing edge of the 
cascade is the same as that at the leading edge. However, in a real engine, the axial velocity can change 
significantly across the cascade due to contraction in the flow path. To assess the effect of streamtube 
contraction, the author derived alternate expressions for the cascade parameters (see Appendix M). 
SetupBFaNS will use these alternate expressions unless Card 5080 in the author file is turned on. 

4.9,5 Equivalent Angle of Attack 

The equivalent angle-of-attack is defined as the angle-of-attack at which an isolated airfoil will have 
the same lift coefficient as the cascade. The equivalent angle-of-attack for an isolated airfoil is determined 
from, 



4.9.6 Upwash 

Turbulent upwash is computed in subroutine “upwash”. To compute upwash, we must transform the 
Reynolds-stress tensor from engine to duct coordinates, and then apply Hanson’s transformation matrix 
(Ref. 4) between duct coordinates and cascade coordinates. Alternately, we can assume isotropic 
turbulence and compute the upwash from the turbulent kinetic energy (this simplified approach is the 
default method in Setup BFaNS, unless Card 5082 is turned on). 

Appendix N shows how Setup BFaNS transforms the Reynolds-stress tensor from the engine 
coordinate system (z,r, 0 ) into the duct coordinate system (i.e., the system that is aligned with the 
streamtubes). The final equations are: 

t 

C\C\ — C\C\ sin 2 (|) + C 3 C 3 cos 2 <j) + 2qc 3 sin (|) cos <j) 

t 

c 2 c 2 =C 2 C 2 

t 

C 3 C 3 = C\C\ cos 2 <() + C 3 C 3 sin 2 (j) - 2 qc 3 sin (|) cos <)) 


C\C 2 — c\C 2 sin <j) + C 7 C 3 cos <() 

C1C3 = ~C\C\ sin <)) cos <|) + C3C3 sin <() cos <|) + qc3 (- cos 2 <)) + sin 2 <()) 

f 

c 2 C 3 = -qc 2 cos (|) + c 2 C 3 sin (|) 

Appendix O shows how Setup BFaNS computes turbulent upwash when Card 5082 is turned on. The 
final equation is: 
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C 2 C 7 ~ C m C n $2n$2n 


where, 


cos — 


I - 5 j cos Y| j's + sin^-| - £ j sin v| j' l sin vp i' s 


-sin] — -£, jcosy'i 


sin £, cosy 5 -cos £, siny^ sinv|/' 5 -cosy£siny' s 


-siny z 


cos] — |cosy'i 


cos £, sin y ' s -sin £, sin y^ cosy 5 sin £, siny' s +cos £, sin \y' L cosy 5 cos y^ cosy 5 


Py = 


The aerodynamic sweep and aerodynamic lean used in the transformation matrix are derived in 
Appendix P. The final equations are: 


Vs = 90 + v|/s - ()> 


tanv| i' L 


tany^ 

-tany^ cos tj) + sin tj) 


4.10 Viscous Quantities 

Setup BFaNS computes the endwall boundary layers at the blade and vane leading edge, in addition 
to the blade-wake thickness at the vane leading edge. Endwall boundary layers are computed with 
subroutines “get bl” and “calc bl”. Blade wakes are computed in “get rotor wake” and 
“calcrotorwake”. 

4.10.1 Endwall Boundary Layers 

The user has the option to manually specify the boundary layer thicknesses in the Setup BFaNS input 
file, or to estimate them according to the procedure described in Appendix Q. 

4.10.2 Rotor Wakes 

Appendix R shows how Setup BFaNS estimates the rotor wakes. 

4.11 Turbulent Length Scale 

Setup BFaNS will default to using the integral length scales that are specified in the flow-field input 
files. However, Cards 5052 or 5054 provide the option to use empirical correlations or the (K, s) values to 
determine the turbulent length scales (see Appendix S). 


5.0 Concluding Remarks 

Pratt & Whitney has developed a Broadband Fan Noise Prediction System (BFaNS) for turbofan 
engines. The noise prediction system consists of two computer programs: SetupBFaNS and BFaNS. 
Setup BFaNS converts user-specified geometry and flow-field information into a BFaNS input file. From 
this input file, BFaNS computes the inlet and aft broadband sound power spectra generated by the fan and 
FEGV. The output file from BFaNS contains the inlet, aft and total sound power spectra from each noise 
source. 

The Broadband Fan Noise Prediction System has been documented in a three-volume report. The 
present volume provides detailed instructions and technical documentation for Setup BFaNS. Volume 2 
provides detailed instructions and technical documentation for BFaNS. Volume 3 provides validation and 
test cases for the prediction system. 
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A.l Variables 


Appendix A. — Symbols 


A 

a 

B 

b 

C or Cj 

C D 

C L 

Cp 

Df 

F 

h 

i 

K 

L 

M 

Pt 

R 

Rl 

R s 

r 

S 

s 

T t 

T,j 

U or Ui 


V 

W or Wi 
x 

y 

z 


area 

speed of sound 
blade count 
airfoil chord 

absolute velocity vector or tensor 
cascade drag coefficient 
cascade lift coefficient 
constant-pressure specific heat 
diffusion factor 
longitudinal spectrum 

cascade gap measured perpendicular to airfoils 
incidence 

turbulent kinetic energy 
length 

Mach number 
total pressure 
gas constant 

Reynolds number based on L 
Reynolds number based on 8 
radius 

streamfunction 
fan-to-FEGV spacing 
total temperature 
generalized 2 nd -order tensor 

blade velocity vector or tensor 
free-stream velocity 
//-component of Reynolds-stress tensor 
vane count 

relative velocity vector or tensor 
horizontal coordinate viewed from rear of engine 
vertical coordinate viewed from rear of engine 
axial coordinate 


A.2 Greek Symbols 

a absolute flow angle; effective angle-of-attack 

P * airfoil metal angle 
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Pm 

b 

Y 

8 

S 2 

8 w 

8 

e 

A 

v 

% 

1 

P 

a 

T 

<t> 

<t>y 

% 

V s 

Vz 

Vs 

Vz 

Vs 

Vz 

Vs 

Q 

Co, Cl, C2 

CO 


relative flow angle 
average flow angle in cascade 

transformation matrix for Cartesian coordinates 

specific-heat ratio 
boundary-layer thickness 
momentum thickness 
wake thickness 
viscous dissipation 
tangential location 
turbulent length scale 
kinematic viscosity 

stagger angle from meridional direction (duct coordinate system) 
stagger angle from tangential direction (cascade coordinate system) 
density 

cascade solidity ( b / r) 
cascade pitch 
meridional flow angle 
energy spectrum tensor 

overall geometric lean angle (engine coordinate system) 

overall geometric sweep angle (engine coordinate system) 
local geometric lean angle (duct coordinate system) 
local geometric sweep angle (duct coordinate system) 
local aerodynamic lean angle (duct coordinate system) 
local aerodynamic sweep angle (duct coordinate system) 
local lean angle (cascade coordinate system) 
local sweep angle (cascade coordinate system) 
fan rotational speed 
coefficients of parabolic fit 
cascade loss coefficient 


A.3 

c 

d 

e 

x 

y 

z 


Subscripts 

cascade coordinates 
duct coordinates 
engine coordinates 
horizontal component 
vertical component 
axial component 
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r radial component 

6 tangential component 

m meridional component; mean cascade 

le leading edge 

te trailing edge 

ID inside diameter 

OD outside diameter 

avg circumferential average 

nose nose-cone leading edge 

high hilite leading edge 

core splitter leading edge 

A.4 Superscripts 

* rotated or translated coordinates 

duct coordinate system 
" cascade coordinate system 
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Appendix B. — Installing the BFaNS System 

The BFaNS system can be downloaded from the following web site: 
http://www.grc.nasa.gov/WWW/5900/5940/code/BFANS/ 

B.l Step Number I: Setting Up the Directory Structure 

1. In your home directory, create a directory called Codes. 

2. Copy the file fans_nasa.tar to your -/Codes directory. 

3. In your -/Codes directory, extract the files using the following command: 

tar -xvf fans_nasa.tar 

4. After executing this command, you will have a ~/Codes/fans_nasa directory. This directory will 
contain the following sub-directories: 

setupbfans: contains source code and make files for Setup BFaNS 
bfans: contains source code and make files for BFaNS 

lib: contains source code and make files for GRAF1C (a freely distributed graphics library created 
by MIT) 

makeflags: contains compiler options for Sun, SGI, and IBM workstations 

5. Appendix K lists the files contained in each directory 

B.2 Step Number II: Setting Up the Graphics Library (GRAFIC) 

1. Go to the ~/Codes/fans_nasa/lib directory. There you will find libraries compiled for an SGI 
workstation (libgrafic.a.Iris), a Sun workstation (libgrafic.a.Sun4), and an IBM workstation 
(libgrafic.a.Ibm6000). You can use one of these existing libraries, or compile the library yourself. If 
you choose to use an existing library, proceed to step no. 2. If you choose to compile the library 
yourself, proceed to step no. 3. 

2. If you are working on an SGI workstation (with iris architecture), type the command: 

cp libgrafic.a.Iris libgrafic.a 

If you are working on a Sun workstation (with svr4 architecture), type the command: 

cp libgrafic.a.Sun4 libgrafic.a 

If you are working on an IBM workstation (with rs6 architecture), type the command: 

cp Iibgrafic.a.lbm6000 libgrafic.a 

For all the above workstations, proceed to Step no. Ill to compile Setup BfaNS. 

3. Go to the ~/Codes/fans_nasa/lib/grafic directory. 

4. If you are working on an SGI workstation (with iris architecture), type the command: 

make -f Makefile. Iris 

If you are working on a Sun workstation (with svr4 architecture), type the command: 

make -f Makefile. Sun4 

If you are working on an IBM workstation (with rs6 architecture), type the command: 

make -f Makefile. Ibm6000 

The above commands will create a file called libgrafic.a in your ~/Codes/fans_nasa/lib directory. If 
you have problems compiling, go to the ~/Codes/fans_nasa/makeflags directory and check the 
compiler options that are set for your workstation (i.e., makeflags.Iris, makeflags.Sun4, or 
makeflags.Ibm6000). 

B.3 Step Number III: Compiling Setup BFaNS 

1 . Go to the ~/Codes/fans_nasa/setup_bfans directory. 

2. If you are working on an SGI workstation (with iris architecture), type the command: 

make -f Makefile. Iris 
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If you are working on a Sun workstation (with svr4 architecture), type the command: 

make -f Makefile. Sun4 

If you are working on an IBM workstation (with rs6 architecture), type the command: 

make -f Makefile. Ibm6000 

The above commands will create an executable file called setupbfans.x in your 
~/Codes/fans_nasa/setup_bfans directory. If you have problems compiling, go to the 
~/Codes/fans_nasa/makeflags directory and check the compiler options that are set for your 
workstation (i.e., makeflags.Iris, makeflags.Sun4, or makeflags.Ibm6000). 

3. If you want to change the compiler optimization (or other options), edit the appropriate makeflags file 
and repeat step no. 2 in this section. 

B.4 Step Number IV: Compiling BFaNS 

1 . Go to the ~/Codes/fans_nasa/bfans directory. 

2. If you are working on an SGI workstation (with iris architecture), type the command: 

make -f Makefile. Iris 

If you are working on a Sun workstation (with svr4 architecture), type the command: 

make -f Makefile. Sun4 

If you are working on an IBM workstation (with rs6 architecture), type the command: 

make -f Makefile. Ibm6000 

The above commands will create an executable file called bfans.x in your ~/Codes/fans_nasa/bfans 
directory. If you have problems compiling, go to the ~/Codes/fans_nasa/makeflags directory and 
check the compiler options that are set for your workstation (i.e., makeflags.Iris, makeflags.Sun4, or 
makeflags.Ibm6000). 

3. If you want to change the compiler optimization (or other options), edit the appropriate makeflags file 
and repeat step no. 2 in this section. 

B.5 Step Number V: Setting Up Aliases 

1. Inside your .cshrc file, create the following aliases: 

alias sbf ~/Codes/fans_nasa/setup_bfans/setup_bfans.x 
alias bf ~/Codes/fans_nasa/bfans/bfans.x 

2. In your home directory, type source .cshrc 

3. These aliases are not required, but they make running much easier. 
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Appendix C. — Running the Test Case 

The BFaNS test case can be downloaded from the following web site: 

http://www.grc.nasa.gov/WWW/5900/5940/code/BFANS/ 

C.l Step Number I: Setting Up the Test Case 

1 . In your home directory, create a directory called bbchallenge. 

2. Copy the file bbchallenge.tar to your -/bbchallenge directory. 

3. In your -/bbchallenge directory, extract the files using the following command: 

tar -xvf bbchallenge.tar 

After executing this command, you will have the following 13 sub-directories inside 
-/bbchallenge: 

bfans_geometry: contains files that describe the flow path, blade, and vane geometry 
bfanscfd: contains files that describe viscous CFD-predicted flow field near blade l.e., blade 
t.e., vane l.e., and vane t.e. for each configuration and operating point 

bfans_data: not used (normally contains files that describe measured flow field near blade l.e., 

blade t.e., vane l.e., and vane t.e. for each configuration and operating point) 

bfans_sline: not used (normally contains files that describe inviscid streamline -predicted flow 

field near blade l.e., blade t.e., vane l.e., and vane t.e. for each configuration and operating point) 

54vr_spdl: contains test case for 54 radial vanes, sideline power 

54vr_spd2: contains test case for 54 radial vanes, cutback power 

54vr_spd3: contains test case for 54 radial vanes, approach power 

26vr_spdl: contains test case for 26 radial vanes, sideline power 

26vr_spd2: contains test case for 26 radial vanes, cutback power 

26vr_spd3: contains test case for 26 radial vanes, approach power 

26vs_spdl: contains test case for 26 radial vanes, sideline power 

26vs_spd2: contains test case for 26 radial vanes, cutback power 

26vs_spd3: contains test case for 26 radial vanes, approach power 

C.2 Step Number II: Running a Test Case for SetupBFaNS 

1. Inside your ~/bbchallenge/54vr_spd3 directory, type the command sbf or 
~/Codes/fans_nasa/setup_bfans/setup_bfans.x. This command will run Setup BFaNS in the current 
directory. 

2. When prompted for the name of the input file, enter setupbfans.input. This file contains the file 
names that describe the geometry and flow field (these files are located in the bfans_geometry and 
bfans_cfd sub-directories under -/bbchallenge), along with some additional information. 

3. A figure will be displayed on the screen showing the blade velocity triangles. Look at the figure and 
confirm that the velocity triangles are correct relative to the blade orientation. Once satisfied, place 
your cursor inside this window, and type the letter x on your keyboard. Typing x will continue 
program execution. 

4. A figure will be displayed on the screen showing the vane velocity triangles. Look at the figure and 
confimi that the velocity triangles are correct relative to the vane orientation. Once satisfied, place 
your cursor inside this window, and type the letter x on your keyboard. Typing x will continue 
program execution. 

5. A figure will be displayed on the screen showing the fan-duct configuration. Look at the figure and 
confimi that the configuration is correct. Once satisfied, place your cursor inside this window, and 
type the letter x on your keyboard. Typing x will continue program execution. 

6. A menu will be displayed to allow you to change the geometry. If you wish to change something, 
type in the number corresponding to your selection and press return. Then provide the requested 
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input. The old and new geometry will now be shown on the screen. Place your cursor inside the figure 
window, and type the letter x on your keyboard. Typing x will continue program execution, allowing 
you to make more changes to the geometry. 

7. Type 99 to finish running Setup BFANS. 

8. You will now have a file named bfans.input in the ~/bbchallenge/54vr_spd3 directory. If you didn’t 
make any changes to the geometry, then bfans.input should be nearly identical to bfans.input.Sun4 
which was created on a Sun workstation at P&W (there might be some very small differences 
depending on the compiler). To confirm that the two files are nearly identical, type the following 
command: 

diff bfans.input bfans.input.Sun4 

This command will show all differences between the two files. 

C.3 Step Number III: Running a Test Case for BFaNS 

1 . Inside your ~/bbchallenge/54vr_spd3 directory, type the command bf or -/Codes/ 
fans_nasa/bfans/bfans.x. This command will run BFaNS inside the current directory (note: to save 
time during code checkout, you might want to edit the bfans.author file to reduce the number of 
frequencies calculated). 

2. When prompted for the name of the input file, enter bfans.input. This file contains the output from 
SetupBFaNS. 

3. When done, BFaNS will display “calculation completed” on the screen. You will now have a file 
named bfans.output in the ~/bbchallenge/54vr_spd3 directory. If you didn’t make any changes to 
the geometry, then bfans.output should be nearly identical to bfans.output.Sun4 which was created 
on a Sun workstation at P&W (there might be some very small differences depending on the 
compiler). To confirm that the two files are nearly identical, type the following command: 

diff bfans.output bfans.output.Sun4 
This command will show all differences between the two files. 
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Appendix D. — Sample Flow-Path Geometry File 


filedescr 

Flowpath: Inside Diameter 

Wall - HUB 


1 

nduct h 

90 

parmname 

Axial 

zhub 

-6.502327 

-6.232568 

-5.962808 

-5.693049 

-5.423289 


-5.153530 

-4.883770 

-4 . 614011 

-4.344251 

-4.074492 


-3.804733 

-3.534973 

-3.265214 

-2.995454 

-2.725695 


-2.455935 

-2 .186176 

-1 . 916416 

-1.646657 

-1.376897 


-1 . 107138 

-0.837379 

-0.567619 

-0.297860 

-0.028100 


0.241659 

0.511419 

0.781178 

1.050938 

1.320697 


1.590456 

1.860216 

2 . 129975 

2.399735 

2 . 669494 


2.939254 

3.209013 

3.478773 

3.748532 

4.018291 


4.288051 

4.557810 

4.827570 

5.097329 

5.367089 


5.636848 

5.906608 

6.176367 

6.446126 

6.715886 


6.985645 

7.255405 

7.525164 

7.794924 

8.064683 


8.334443 

8.604202 

8.873961 

9.143721 

9.413480 


9.683240 

9.952999 

10.222759 

10.492518 

10.762278 


11.032037 

11.301796 

11.571556 

11 . 841315 

12 .111075 


12.380834 

12.650594 

12.920353 

13.190113 

13.459872 


13.729631 

13.999391 

14.269150 

14.538910 

14.808669 


15.078429 

15.348188 

15.617948 

15.887707 

16.157467 


16.427226 

16.696985 

16.966745 

17.236504 

17.506264 

parmname 

Radius | 

rhub 

2.250805 

2.250757 

2.250756 

2.250848 

2.250945 


2.250903 

2.250703 

2.251318 

2.250003 

2.268878 


2.297416 

2.373692 

2.525183 

2 . 674208 

2.791127 


2.900357 

3.002891 

3.099274 

3.189671 

3.276292 


3.358367 

3.435867 

3.509917 

3.581078 

3.650682 


3.720369 

3.794309 

3.878063 

3.969007 

4.067281 


4 . 170762 

4.275789 

4.379734 

4.479830 

4.572974 


4.656779 

4.731721 

4.799961 

4.862569 

4 . 921029 


4.976265 

5.029055 

5.079441 

5.127756 

5.173022 


5.214311 

5.251548 

5.287234 

5.322437 

5.356940 


5.391048 

5.429903 

5.478404 

5.534245 

5.592687 


5.649140 

5.699345 

5.740850 

5.772924 

5.795676 


5.810992 

5.821017 

5.827340 

5.831110 

5.833467 


5.835088 

5.836256 

5.836969 

5.836978 

5.836439 


5.836241 

5.837039 

5.839078 

5.843604 

5.850944 


5.859892 

5.875535 

5.892634 

5.916299 

5.941113 


5.971357 

6.001601 

6.037731 

6.073982 

6.113309 


6.155001 

6.196693 

6.242894 

6.289515 

6.336137 

blank 

1 

filedesc 

Primary Splitter: 

Inside Diameter Wall - 

UPPER 

1 

nsplt u 

50 

parmname 

Axial 

zsup 

6.756294 

6.975681 

7.195069 

7.414456 

7.633843 


7.853230 

8.072617 

8.292004 

8.511391 

8.730778 


8.950166 

9.169553 

9.388940 

9.608327 

9.827714 


10.047101 

10.266488 

10.485875 

10.705263 

10.924650 


11 . 144037 

11.363424 

11.582811 

11.802198 

12.021585 


12.240973 

12.460360 

12 . 679747 

12 . 899134 

13.118521 


13.337908 

13.557295 

13.776682 

13.996070 

14.215457 


14.434844 

14 . 654231 

14.873618 

15.093005 

15.312392 


15.531780 

15.751167 

15.970554 

16.189941 

16.409328 


16.628715 

16.848102 

17.067489 

17.286877 

17.506264 

parmname 

Radius | 

rsup 

5.362028 

5.389747 

5.420434 

5.457386 

5.500269 


5.546825 

5.594395 

5.640565 

5.683004 

5.719951 


5.750893 

5.775476 

5.793916 

5.807358 

5.816899 


5.823565 

5.828091 

5.831034 

5.833065 

5.834527 


5.835635 

5.836469 

5.836996 

5.836986 

5.836640 


5.836273 

5.836223 

5.837216 

5.838722 

5.842403 


5.846898 

5.854175 

5.862536 

5.875341 

5.888147 


5.907170 

5.926415 

5.948395 

5.972991 

5.997588 


6.026151 

6.055633 

6.085115 

6.118328 

6.152235 


6.186142 

6.222389 

6.260305 

6.298221 

6.336137 

blank 

1 

filedesc 

Primary Splitter: 

Outside 

Diameter Wall 

- LOWER 

1 


N ASA/CR— 20 1 0-2 1 6898/VOL 1 


37 




nsplt 1 

50 

parmname 

Axial 

zslo 

6.756294 

6.975681 

7.195069 

7.414456 

7.633843 


7.853230 

8.072617 

8.292004 

8.511391 

8.730778 


8.950166 

9.169553 

9.388940 

9.608327 

9.827714 


10.047101 

10.266488 

10.485875 

10.705263 

10.924650 


11 . 144037 

11.363424 

11.582811 

11 . 802198 

12.021585 


12.240973 

12.460360 

12 . 679747 

12.899134 

13.118521 


13.337908 

13.557295 

13.776682 

13.996070 

14.215457 


14.434844 

14 . 654231 

14.873618 

15.093005 

15.312392 


15.531780 

15.751167 

15.970554 

16.189941 

16.409328 


16.628715 

16.848102 

17.067489 

17.286877 

17.506264 

parmname 

Radius | 

rslo 

5.362028 

5.389747 

5.420434 

5.457386 

5.500269 


5.546825 

5.594395 

5.640565 

5.683004 

5.719951 


5.750893 

5.775476 

5.793916 

5.807358 

5.816899 


5.823565 

5.828091 

5.831034 

5.833065 

5.834527 


5.835635 

5.836469 

5.836996 

5.836986 

5.836640 


5.836273 

5.836223 

5.837216 

5.838722 

5.842403 


5.846898 

5.854175 

5.862536 

5.875341 

5.888147 


5.907170 

5.926415 

5.948395 

5.972991 

5.997588 


6.026151 

6.055633 

6.085115 

6.118328 

6.152235 


6.186142 

6.222389 

6.260305 

6.298221 

6.336137 

blank 

1 

filedesc 

Flowpath: Outside Diameter 

Wall - TIP 


1 

nduct t 

35 

parmname 

Axial 

ztip 

-6.502327 

-6.232568 

-5.962808 

-5.693049 

-5.423289 


-5.153530 

-4.883770 

-4 . 614011 

-4.344251 

-4.074492 


-3.804733 

-3.534973 

-3.265214 

-2.995454 

-2.725695 


-2.455935 

-2 .186176 

-1 . 916416 

-1.646657 

-1.376897 


-1 . 107138 

-0.837379 

-0.567619 

-0.297860 

-0.028100 


0.241659 

0.511419 

0.781178 

1.050938 

1.320697 


1.590456 

1.860216 

2 . 129975 

2.399735 

2 . 669494 


2.939254 

3.209013 

3.478773 

3.748532 

4.018291 


4.288051 

4.557810 

4.827570 

5.097329 

5.367089 


5.636848 

5.906608 

6.176367 

6.446126 

6.715886 


6.985645 

7.255405 

7.525164 

7.794924 

8.064683 


8.334443 

8.604202 

8.873961 

9.143721 

9.413480 


9.683240 

9.952999 

10.222759 

10.492518 

10.762278 


11.032037 

11.301796 

11.571556 

11 . 841315 

12 .111075 


12.380834 

12.650594 

12.920353 

13.190113 

13.459872 


13.729631 

13.999391 

14.269150 

14.538910 

14.808669 


15.078429 

15.348188 

15.617948 

15.887707 

16.157467 


16.427226 

16.696985 

16.966745 

17.236504 

17.506264 

parmname 

Radius | 

rtip 

11.010941 

11.012493 

11.013998 

11.015268 

11.016401 


11.017396 

11.018124 

11.018691 

11.019068 

11.019232 


11.019168 

11.018869 

11.018327 

11.017558 

11.016528 


11.015241 

11.013677 

11.011839 

11.011221 

11.006656 


10.995469 

10.981015 

10.966768 

10.953494 

10.940141 


10.926359 

10.912204 

10.898097 

10.884698 

10.873096 


10.865967 

10.865466 

10.870844 

10.879611 

10.891271 


10.905508 

10.922375 

10.941645 

10.962651 

10.984550 


11.006749 

11.028477 

11.048728 

11.067650 

11.085357 


11 . 101904 

11.117003 

11.130196 

11 . 141621 

11.150951 


11.157653 

11.160125 

11.156557 

11 . 145514 

11 . 127017 


11 . 104327 

11.080913 

11.059781 

11.042190 

11.028078 


11.017560 

11.010750 

11.006802 

11.004715 

11.003883 


11.003863 

11.003966 

11.003947 

11.003903 

11.003971 


11.003982 

11.003810 

11.003820 

11.006583 

11.013209 


11.022273 

11.036808 

11.051996 

11.069672 

11.087250 


11 . 104318 

11 . 121385 

11 . 134995 

11 . 148494 

11.158052 


11 .164576 

11 . 171100 

11.168430 

11.164928 

11 .161425 


Note: the shaded column is not part of the input file 
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D.l FORTRAN Code Used to Read File 


read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit, 
read (unit. 


' (a50) 

') 

filedescr 


*) 


nduct h 


' (a50) 

') 

parmname 


*) 


(zhub ( i) , i=l , nduct 

_h) 

' (a50) 

') 

parmname 


*) 


(rhub ( i) , i=l , nduct 

_h) 

' (al) ' 

) 

blank 


' (a50) 

') 

filedescr 


*) 


nsplt u 


' (a50) 

') 

parmname 


*) 


(zsup (i) , i=l , nsplt 

_u) 

' (a50) 

') 

parmname 


*) 


(rsup ( i) , i=l , nsplt 

_u) 

' (al) ' 

) 

blank 


' (a50) 

') 

filedescr 


*) 


nsplt 1 


' (a50) 

') 

parmname 


*) 


(zslo(i) , : — 1 , nsplt 

_1> 

' (a50) 

') 

parmname 


*) 


(rslo(i) , i=l, nsplt 

_1) 

' (al) ' 

) 

blank 


' (a50) 

') 

filedescr 


*) 


nduct t 


' (a50) 

') 

parmname 


*) 


(ztip ( i) , i=l , nduct 

_t) 

' (a50) 

') 

parmname 


*) 


(rtip ( i) , i=l , nduct 

_t) 


hub surface description 
number of points along hub 
hub axial coordinate name 
hub axial coordinates 
hub radial coordinate name 
hub radial coordinates 
blank line 

upper splitter surface description 
number of points along upper splitter 
upper splitter axial coordinate name 
upper splitter axial coordinates 
upper splitter radial coordinate name 
upper splitter radial coordinates 
blank line 

lower splitter surface description 
number of points along lower splitter 
lower splitter axial coordinate name 
lower splitter axial coordinates 
lower splitter radial coordinate name 
lower splitter radial coordinates 
blank line 

shroud surface description 
number of points along shroud 
shroud axial coordinate name 
shroud axial coordinates 
shroud radial coordinate name 
hub radial coordinates 
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Appendix E. — Sample Blade Geometry File 


filedescr 

NASA SOURCE DIAGNOSTICS FAN 

ROTOR AIRFOIL 

DATA 

1 

numrle 

51 

parmname 

XLE 

xle 

-3.300717 

-3.303915 

-3.308069 

-3.313466 

-3.320478 


-3.329237 

-3.340175 

-3.353832 

-3.370878 

-3.392150 


-3.418686 

-3.451780 

-3.493053 

-3.544538 

-3.608829 


-3.689264 

-3.790090 

-3.916614 

-4.075167 

-4.273249 


-4.520267 

-4.828773 

-5.215665 

-5.700811 

-6.305202 


-7.061386 

-7.823350 

-8.427943 

-8.918634 

-9.317897 


-9.638566 

-9.894466 

-10.098216 

-10.260604 

-10.390247 


-10.493849 

-10.576728 

-10.643064 

-10.696157 

-10.738645 


-10.772639 

-10.799839 

-10.821598 

-10.839005 

-10.852930 


-10.864068 

-10.893050 

-10.872978 

-10.880105 

-10.885588 

-10.889806 

parmname 

YLE | 

yle 

0.056513 

0.058546 

0.061189 

0.064622 

0.069084 


0.074656 

0.081611 

0.090283 

0.101076 

0.114469 


0.131012 

0.151285 

0.175824 

0.204982 

0.238555 


0.275329 

0.313901 

0.353320 

0.395976 

0.448298 


0.516396 

0.600231 

0.691755 

0.792850 

0.932052 


1.099041 

1.227651 

1.360980 

1.423135 

1.429222 


1.423645 

1.421848 

1.426293 

1.433713 

1.441308 


1.448014 

1.453255 

1.457092 

1.459912 

1.462011 


1.463594 

1.464802 

1.465734 

1.466459 

1.467026 


1.467472 

1.468602 

1.467826 

1.468104 

1.468316 

1.468480 

parmname 

ZLE | 

zle 

-1.297397 

-1.297990 

-1.298759 

-1.299756 

-1.301048 


-1.302654 

-1.304649 

-1.307122 

-1.310177 

-1.313935 


-1.318526 

-1.324082 

-1.330689 

-1.338260 

-1.346318 


-1.353738 

-1.359192 

-1.361348 

-1.360299 

-1.357423 


-1.351905 

-1.338306 

-1.312679 

-1.274927 

-1.226397 


-1.165348 

-1 .119212 

-1 . 100427 

-1.086893 

-1.073364 


-1.064533 

-1.060213 

-1.057601 

-1.055210 

-1.052501 


-1.049508 

-1.046523 

-1.043761 

-1.041353 

-1.039327 


-1.037654 

-1.036288 

-1.035179 

-1.034283 

-1.033560 


-1.032978 

-1.031453 

-1.032511 

-1.032136 

-1.031847 

-1.031624 

parmname 

BE TATI | 

betale 

66.841801 

66.815885 

66.783103 

66.742593 

66.683989 


66.614494 

66.518115 

66.399710 

66.240933 

66.033458 


65.750624 

65.367511 

64.835773 

64 .115349 

63.168008 


61 . 994470 

60.685418 

59.429298 

58.195041 

56.664610 


54.802715 

53.044276 

51.019784 

48.935041 

46.509141 


43.978135 

41 . 697607 

39.574998 

38.087785 

36.851772 


35.808748 

35.100205 

34.436097 

33.761993 

33.149026 


32.648950 

32.246985 

31.915624 

31.643311 

31.425093 


31.245763 

31 . 101723 

30.986595 

30.895777 

30.819186 


30.759661 

30.608167 

30.714022 

30.676131 

30.646688 

30.624129 

parmname 

XTE | 

xte 

-4.088949 

-4.092535 

-4.097195 

-4.103248 

-4.111110 


-4 . 120854 

-4 . 132927 

-4 . 147878 

-4.166385 

-4 .189279 


-4.217576 

-4.252514 

-4.295610 

-4.348718 

-4.414104 


-4.494511 

-4.593260 

-4.714419 

-4.862988 

-5.045069 


-5.268126 

-5.541716 

-5.878458 

-6.292905 

-6.801577 


-7.431083 

-8.067842 

-8.579360 

-8.995760 

-9.339889 


-9.620321 

-9.846050 

-10.027071 

-10.172226 

-10.288950 


-10.382961 

-10.458827 

-10.520140 

-10.569674 

-10.609668 


-10.641939 

-10.667974 

-10.688969 

-10.705896 

-10.719543 


-10.730546 

-10.759545 

-10.739413 

-10.746562 

-10.752061 

-10.756291 

parmname 

YTE | 
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yte 

0.595555 

0.593313 

0.590397 

0.586604 

0.581671 


0.575546 

0.567943 

0.558513 

0.546824 

0.532360 


0.514516 

0.492602 

0.465780 

0.433002 

0.392969 


0.344221 

0.285327 

0.214571 

0.129390 

0.026424 


-0.097081 

-0.237919 

-0.392177 

-0.576410 

-0.812905 


-1.075679 

-1.269325 

-1.430582 

-1.533526 

-1.565878 


-1.572818 

-1.579427 

-1.590009 

-1.603047 

-1.615155 


-1.625258 

-1 . 632740 

-1 . 637674 

-1.640837 

-1 . 642914 


-1 . 644310 

-1 . 645274 

-1 . 645954 

-1 . 646446 

-1.646806 


-1 . 647074 
-1.647700 

-1 . 647280 

-1 . 647433 

-1.647550 

-1.647635 

parmname 

ZTE | 

zte 

1.490532 

1.491784 

1.493412 

1.495524 

1.498265 


1.501657 

1.505851 

1.511032 

1.517420 

1.525280 


1.534918 

1.546681 

1.560933 

1.578008 

1.598115 


1 . 621252 

1 . 647092 

1 . 675143 

1.705382 

1.738630 


1.773224 

1.795254 

1.785352 

1.761071 

1.718796 


1.631256 

1.531697 

1.440159 

1.362147 

1.302295 


1.252064 

1.208955 

1 . 174382 

1 . 147466 

1 . 126418 


1.109854 

1.096681 

1.086109 

1.077594 

1.070730 


1.065192 

1.060726 

1.057124 

1.054220 

1.051879 


1.049991 

1.045016 

1.048470 

1.047244 

1.046300 

1.045575 

parmname 

BE TAT 2 | 

betate 

123.481562 

123.359068 

123.200406 

122 . 994211 

122.723927 


122.390499 

121 . 977237 

121.463885 

120.829013 

120.048534 


119.084911 

117 . 907573 

116.482608 

114.764199 

112.735775 


110.376564 

107.716961 

104.740661 

101.130447 

96.493614 


90.782852 

84.861986 

79.986424 

74.797613 

68.264310 


59.838997 

53.438499 

49.173076 

45.258556 

43.025379 


42.047168 

41.267618 

40.301256 

39.256656 

38.427494 


37 . 915711 

37.615648 

37.431940 

37.300851 

37.201613 


37 . 121749 

37.054726 

37.000575 

36.955608 

36.922670 


36.894472 

36.812149 

36.866245 

36.849647 

36.828629 

36.820108 


Note: the shaded column is not part of the input file 


E.l FORTRAN Code Used to Read File 


read (unit, ' (a50) ' ) 

f iledescr 

1 

file 

description 

read (unit, *) 

numrle 

1 

number of points defini 

read(unit, ' (a50) ') 

parmname 

1 

CD 

i — 1 

x-coordinate name 

read (unit, *) 

(xle (i) , i=l, numrle) 

1 

CD 

i — 1 

x-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

<D 

i — 1 

y-coordinate name 

read (unit, *) 

(yle (i) , i=l, numrle) 

1 

CD 

i — 1 

y-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

CD 

i — 1 

z-coordinate name 

read (unit, *) 

(zle (i) , i=l, numrle) 

1 

1 . e . 

z-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

Q) 

i — i 

metal angle name 

read (unit, *) 

(betale (i) , i=l, numrle) 

1 

CD 

i — 1 

metal angle 

read (unit, ' (a50) ' ) 

parmname 

1 

t . e . 

x-coordinate name 

read (unit, *) 

(xte (i) , i— 1, numrle) 

1 

t . e . 

x-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

t . e . 

y-coordinate name 

read (unit, *) 

(yte (i) , i=l, numrle) 

1 

t . e . 

y-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

t . e . 

z-coordinate name 

read (unit, *) 

(zte (i) , i=l, numrle) 

| 

t . e . 

z-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

t . e . 

metal angle name 

read (unit, *) 

(betate (i) , i=l, numrle) 

1 

t . e . 

metal angle 
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Appendix F. — Sample Vane Geometry File 


filedescr 

NASA SOURCE DIAGNOSTICS VANE 

STATOR AIRFOIL 

DATA 

1 

numrle 

51 

parmname 

XLE 

xle 

-5.416727 

-5.419984 

-5.424215 

-5.429717 

-5.436869 


-5.445656 

-5.456452 

-5.469716 

-5.486012 

-5.506033 


-5.530631 

-5.560848 

-5.597970 

-5.643573 

-5.699591 


-5.768402 

-5.852918 

-5.956717 

-6.084190 

-6.240739 


-6.433002 

-6.669145 

-6.959180 

-7.315404 

-7.752908 


-8.290226 

-8.827341 

-9.264373 

-9.619956 

-9.909304 


-10.144755 

-10.336369 

-10.492330 

-10.619282 

-10.722619 


-10.806724 

-10.875172 

-10.930878 

-10.976212 

-11.013106 


-11.043136 

-11.067576 

-11.087468 

-11.103658 

-11.116835 


-11 . 127559 
-11.156293 

-11 . 136287 

-11 . 143391 

-11.148856 

-11.153060 

parmname 

YLE | 

yle 

-0.406363 

-0.406194 

-0.405973 

-0.405686 

-0.405312 


-0.404852 

-0.404287 

-0.403591 

-0.402734 

-0.401682 


-0.400392 

-0.398819 

-0.396905 

-0.394579 

-0.391767 


-0.388400 

-0.384440 

-0.379891 

-0.374796 

-0.369099 


-0.362613 

-0.354976 

-0.345888 

-0.335125 

-0.322678 


-0.308797 

-0.296720 

-0.288643 

-0.284890 

-0.283530 


-0.283985 

-0.285139 

-0.285764 

-0.285726 

-0.285371 


-0.285092 

-0.285040 

-0.285216 

-0.285526 

-0.285869 


-0.286190 

-0.286472 

-0.286714 

-0.286918 

-0.287089 


-0.287229 

-0.287618 

-0.287347 

-0.287443 

-0.287516 

-0.287574 

parmname 

ZLE | 

zle 

7.267998 

7.268016 

7.268040 

7.268072 

7.268112 


7.268163 

7.268225 

7.268301 

7.268394 

7.268508 


7.268649 

7.268822 

7.269035 

7.269298 

7.269619 


7.270013 

7.270496 

7.271083 

7.271793 

7.272658 


7.273715 

7.275023 

7.276612 

7.278489 

7.280753 


7.283655 

7.286515 

7.288820 

7.290805 

7.292380 


7.293630 

7.294653 

7.295467 

7.296127 

7.296694 


7.297163 

7.297535 

7.297823 

7.298047 

7.298222 


7.298362 

7.298473 

7.298564 

7.298636 

7.298696 


7.298744 

7.298871 

7.298783 

7.298814 

7.298839 

7.298857 

parmname 

BE TATI | 

betale 

131.281769 

131.272356 

131.246212 

131.223061 

131.190830 


131 . 152542 

131.100057 

131.044341 

130.971196 

130.874959 


130.761550 

130.623424 

130.456048 

130.250116 

129.982424 


129.669266 

129.283355 

128.854058 

128.385733 

127 . 861015 


127.242161 

126.511595 

125.711721 

124.809808 

123.857544 


123.024853 

122 . 618754 

122.587236 

122.795412 

123.133131 


123.582709 

124.028165 

124.490860 

124.887123 

125.190549 


125.514202 

125.867713 

126.185061 

126.509030 

126.772852 


127.005851 

127 .182839 

127.330880 

127.471321 

127.568032 


127 . 651582 
127.867114 

127.713520 

127.773353 

127 . 810731 

127 . 848052 

parmname 

XTE | 

xte 

-5.731794 

-5.735029 

-5.739233 

-5.744697 

-5.751801 


-5.760493 

-5.771126 

-5.784135 

-5.800049 

-5.819518 


-5.843336 

-5.872474 

-5.908123 

-5.951734 

-6.005087 


-6.070358 

-6.150211 

-6.247900 

-6.367413 

-6.513622 


-6.692492 

-6.911319 

-7 . 179029 

-7.506539 

-7 . 907212 


-8.397385 

-8.887561 

-9.288233 

-9.615743 

-9.883452 


-10.102280 

-10.281151 

-10.427359 

-10.546871 

-10.644561 


-10.724414 

-10.789686 

-10.843039 

-10.886650 

-10.922299 


-10.951438 

-10.975257 

-10.994726 

-11.010640 

-11.023649 


-11.034283 

-11.062981 

-11.042974 

-11.050076 

-11.055543 

-11.059748 

parmname 

YTE 1 
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yte 

0.000471 

0.000471 

0.000471 

0.000472 

0.000471 


0.000470 

0.000470 

0.000469 

0.000469 

0.000468 


0.000469 

0.000470 

0.000471 

0.000473 

0.000473 


0.000470 

0.000468 

0.000469 

0.000473 

0.000473 


0.000472 

0.000471 

0.000468 

0.000475 

0.000471 


0.000464 

0.000464 

0.000468 

0.000471 

0.000460 


0.000467 

0.000467 

0.000462 

0.000457 

0.000460 


0.000466 

0.000469 

0.000470 

0.000467 

0.000464 


0.000463 

0.000463 

0.000463 

0.000463 

0.000466 


0.000465 

0.000468 

0.000467 

0.000467 

0.000467 

0.000468 

parmname 

ZTE | 

zte 

8.809691 

8.809727 

8.809773 

8.809833 

8.809911 


8.810006 

8.810122 

8.810264 

8.810437 

8.810648 


8.810906 

8.811219 

8.811600 

8.812058 

8.812607 


8.813257 

8.814007 

8.814852 

8.815793 

8.816864 


8.818124 

8.819621 

8.821311 

8.823187 

8.825278 


8.827568 

8.829410 

8.830547 

8.831099 

8.831244 


8.831086 

8.830804 

8.830573 

8.830452 

8.830420 


8.830431 

8.830447 

8.830440 

8.830409 

8.830365 


8.830317 

8.830269 

8.830224 

8.830184 

8.830150 


8.830121 

8.830036 

8.830096 

8.830075 

8.830059 

8.830048 

parmname 

BE TAT 2 | 

betate 

86.977677 

86.973932 

86.976088 

86.973618 

86.968883 


86.975419 

86.969165 

86.970769 

86.963695 

86.964786 


86.966653 

86.968617 

86.977814 

86.980755 

86.998967 


86.987708 

86.995244 

86.996859 

87.022108 

87.024978 


87.028719 

87.050052 

87.087533 

87 .184788 

87.230916 


87.323123 

87.347619 

87.449030 

87.501262 

87.537642 


87.590556 

87.606019 

87.640125 

87.669546 

87.656125 


87 . 649347 

87 . 658172 

87.670551 

87.678676 

87 . 676479 


87.700663 

87.709592 

87.709429 

87.701384 

87.694857 


87 . 696417 
87.682160 

87.686811 

87.693660 

87 . 671975 

87 . 677902 


Note: the shaded column is not part of the input file 


F.l FORTRAN Code Used to Read File 


read (unit, ' (a50) ' ) 

f iledescr 

1 

file 

description 

read (unit, *) 

numrle 

1 

number of points defini 

read(unit, ' (a50) ') 

parmname 

1 

CD 

i — 1 

x-coordinate name 

read (unit, *) 

(xle (i) , i=l, numrle) 

1 

CD 

i — 1 

x-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

<D 

i — 1 

y-coordinate name 

read (unit, *) 

(yle (i) , i=l, numrle) 

1 

CD 

i — 1 

y-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

CD 

i — 1 

z-coordinate name 

read (unit, *) 

(zle (i) , i=l, numrle) 

1 

1 . e . 

z-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

Q) 

i — i 

metal angle name 

read (unit, *) 

(betale (i) , i=l, numrle) 

1 

CD 

i — 1 

metal angle 

read (unit, ' (a50) ' ) 

parmname 

1 

t . e . 

x-coordinate name 

read (unit, *) 

(xte (i) , i=l, numrle) 

1 

t . e . 

x-coordinates 

read (unit, ' (a50) ' ) 

string 

1 

t . e . 

y-coordinate name 

read (unit, *) 

(yte (i) , i=l, numrle) 

1 

t . e . 

y-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

t . e . 

z-coordinate name 

read (unit, *) 

(zte (i) , i=l, numrle) 

| 

t . e . 

z-coordinates 

read (unit, ' (a50) ' ) 

parmname 

1 

t . e . 

metal angle name 

read (unit, *) 

(betate (i) , i=l, numrle) 

1 

t . e . 

metal angle 
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Appendix G. — FORTRAN Code Used to Read Flow-Field File 
(Viscous CFD and Experimental Data) 

read (unitnum, * ) iwaketyp 

read (unitnum, * ) omega, xprtip, rhub, rtip, npassage, phipr 

read (unitnum, * ) numr,numt 

read (unitnum, * ) ( z ( i ) , i=l , numr ) 

read (unitnum, * ) ((x(i,j),i=l, numr) , j=l , numt) 

read (unitnum, * ) ((y(i,j),i=l, numr) , j=l , numt) 

read (unitnum,*) icz,icu,icr,itke,ieps,icpzcpz, icpucpu, 

& icprcpr, icpzcpu, icpzcpr, icpucpr , iscale 

if (icz .eq. 1) then 

read (unitnum, *) ( (cz (i, j ) , i=l, numr) , j = l, numt) 

endif 

if (icu .eq. 1) then 

read (unitnum, *) ( (cu ( i, j ) , i=l , numr) , j=l, numt) 

endif 

if (icr .eq. 1) then 

read (unitnum, *) ((cr(i,j),i=l, numr) , j=l , numt) 

endif 

if (itke .eq. 1) then 

read (unitnum, *) ( (tke (i, j ) , i=l, numr) , j =1 , numt) 

endif 

if (ieps .eq. 1) then 

read (unitnum, *) ( (eps (i, j ) , i=l, numr) , j =1 , numt) 

endif 

if (icpzcpz .eq. 1) then 

read (unitnum, *) ( (cpzcpz ( i, j ) , i=l, numr) , j=l, numt) 

endif 

if (icpucpu .eq. 1) then 

read (unitnum, *) ( (cpucpu ( i, j ) , i=l, numr) , j=l, numt) 

endif 

if (icprcpr .eq. 1) then 

read (unitnum, *) ( (cprcpr ( i, j ) , i=l, numr) , j=l, numt) 

endif 

if (icpzcpu .eq. 1) then 

read (unitnum, * ) ( (cpzcpu ( i, j ) , i=l , numr ) , j =1 , numt) 

endif 

if (icpzcpr .eq. 1) then 

read (unitnum, *) ( (cpzcpr ( i, j ) , i=l, numr) , j=l, numt) 

endif 

if (icpucpr .eq. 1) then 

read (unitnum, *) ( (cpucpr ( i, j ) , i=l, numr) , j=l, numt) 

endif 

if (iscale .eq. 1) then 

read (unitnum, *) (scale (i) , i=l,numr) 

endif 

read (unitnum,*) ztip2 , rhub2 , rtip2 

read (unitnum,*) numr2 

read (unitnum,*) ( r2 ( i ) , i=l , numr2 ) 

read (unitnum,*) ( z2 ( i ) , i=l , numr2 ) 

read (unitnum,*) (pt2 ( i ) , i=l , numr2 ) 

read (unitnum,*) ( tt2 ( i ) , i=l , numr2 ) 
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Appendix H. — FORTRAN Code Used to Read Flow-Field File 
(Inviscid Streamline Code) 


read (unitnum, ' (a50) 
read (unitnum,*) 
read (unitnum, ' (a50) 
read (unitnum,*) 
read (unitnum, ' (a50) 
read (unitnum,*) 
read (unitnum, ' (a50) 
read (unitnum,*) 
do i = 1 , numr 

r (i) = diam(i) /2 . 
z(i) = z(i) - zblade/in 


enddo 




numt 

- 1 



rhub 

= r (1) 



rtip 

= r (numr) 



omega 

= (u (1) *12 . ) /r (1) 


read 

(unitnum, ' (a50) 

') 

parmname 

read 

(unitnum, *) 


( (cz (i, j ) , i=l , numr) , j=l , numt) 

read 

(unitnum, ' (a50) 

') 

parmname 

read 

(unitnum, *) 


( (ct (i, j ) , i=l, numr) , j=l , numt) 

read 

(unitnum, ' (a50) 

') 

parmname 

read 

(unitnum, *) 


( (cr (i, j ) , i=l, numr) , j=l , numt) 

read 

(unitnum, ' (a50) 

') 

parmname 

read 

(unitnum, * ) 


(loss (i) , i=l , numr) 

read 

(unitnum, ' (a50) 

') 

parmname 

read 

(unitnum, *) 


(gapchord (i) , i=l, numr) 

numr2 

= numrl 



rhub2 

= rhubl/in 



rtip2 

= rtipl/in 



read 

(unitnum, ' (a50) 

') 

parmname 

read 

(unitnum, *) 


(pt2 (i) , 1=1, numr 2) 

read 

(unitnum, ' (a50) 

') 

parmname 

read 

(unitnum, *) 


(tt2 (i) , i=l, numr 2 ) 

do i 

= l,numr2 




r2 (i) = rl (i) /in 
enddo 


' ) filedescr 
numr 

' ) parmname 

(diam ( i ) , i=l , numr ) 
' ) parmname 

(z (i) , i=l, numr) 

' ) parmname 

(u ( i ) , i=l , numr) 
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Appendix I. — Sample Setup BFaNS Input File 


case_54vr_ap . data 

1.0450 

8.8300 

-6.5023 

-4 .0074 

6.7563 

17.5063 

rotor_54vr_ap. data 
22 

-1 

stator_54vr_ap . data 
54 

wake_54vr_ap_rotorle . data 
2 

wake_54vr_ap_rotorte . data 
2 

wake_54vr_ap_statorle . data 
2 

wake_54vr_ap_statorte . data 
2 

0.06 

0.09 

0.13 

0.64 

54-Vane Radial Stator 
Approach Condition (7,808 rpm) 


10. DUCT GEOMETRY FILE 

20. AXIAL COORDINATE OF BLADE TRAILING EDGE TIP (IN) 
30. AXIAL COORDINATE OF VANE TRAILING EDGE TIP (IN) 
40. AXIAL COORDINATE OF HILITE (IN) 

50. AXIAL COORDINATE OF NOSE CONE (IN) 

60. AXIAL COORDINATE OF CORE INLET (IN) 

70. AXIAL COORDINATE OF DUCT EXIT (IN) 

80. ROTOR GEOMETRY FILE 
90. BLADE NUMBER 

100. DIRECTION OF ROTATION FROM REAR (1: CCW, -1: CW) 
110. STATOR GEOMETRY FILE 
120. VANE NUMBER 

130. ROTOR UPSTREAM FLOWFIELD FILE 

140. PRECEDING FILE TYPE (1=SLINE, 2=CFD, 3=EXPER) 
150. ROTOR DOWNSTREAM FLOWFIELD FILE 
160. PRECEDING FILE TYPE (1=SLINE, 2=CFD, 3=EXPER) 
170. STATOR UPSTREAM FLOWFIELD FILE 
180. PRECEDING FILE TYPE (1=SLINE, 2=CFD, 3=EXPER) 
190. STATOR DOWNSTREAM FLOWFIELD FILE 
200. PRECEDING FILE TYPE (1=SLINE, 2=CFD, 3=EXPER) 
210. B.L. THKNESS AT BLADE HUB (IN) (0 FOR AUTOCALC) 

220. B.L. THKNESS AT BLADE TIP (IN) (0 FOR AUTOCALC) 

230. B.L. THKNESS AT VANE HUB (IN) (0 FOR AUTOCALC) 

240. B.L. THKNESS AT VANE TIP (IN) (0 FOR AUTOCALC) 

250. RUN TITLE #1 
260. RUN TITLE #2 


Notes : 


1. Numeric input is free format. 

2. Character input is formatted as A30, and must not be in quotes. 

3. The first 30 columns are reserved for user input. 

4. Comments may be added in columns 31 and greater. 

5. The blade and vane number must correspond to the blade and vanes numbers 
used to generate the flow field files 
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Appendix J. — Sample Setup BFaNS Author File 


c 

c author file for setup_bfans 

e 

c section #1: read format is ( al , i5, a2 0, a54 ) 

c 

c card description : a54 

c 


-5000 input file name 
-5002 geometry dir 
-5004 streamline dir 
-5006 cfd dir 
-5008 data dir 
-5010 source dir 
-5012 through 5048 


setupbfans . input 
/bfans_geometry/ 

/bfans_sline/ 

/bfans_cfd/ 

/bfans_data/ 

~ /Codes /fans_nasa/ setup_bfans/ 
reserved for future development 


section #2: read format is (al , i5, a4 5, 2i4 , 2f 8 , a5) 




card description : i4 i4 

f 8 

f 8 

a5 


-5050 create debug files 

-5052 do not use measured scale at station 
-5054 do not use measured scale at station 
-5056 over-ride scale * epsilon / tke**1.5 
-1 defaults to 1.00 
-5058 tke fudge factor 
-1 defaults to 1.00 
-5060 epsilon fudge factor: 

-1 defaults to (tke fudge factor) **1.5 
-5062 over-ride scale cal. factor at blade 
-1 rmoda = cal. factor (default = 0.62) 

-5064 over-ride scale cal. factor at blade 
-1 rmoda = cal. factor (default = 0.62) 

-5066 over-ride scale cal. factor at vane hub : 0.62 

-1 rmoda = cal. factor (default =0.62) : 

-5068 over-ride scale cal. factor at vane tip : 0.62 

-1 rmoda = cal. factor (default =0.62) : 

-5070 over-ride scale cal. factor for blade wake : 0.68 

-1 rmoda = cal. factor (default =0.68) : 

-5072 over-ride order of wake curve fit : 5 

-1 default = 5 : 

-5074 over-ride tol . for unwrapping blade wake : 0.90 

-1 rmoda = tolerance (def =0.9) : 

-5076 write wake lean to output file : 

-5078 account for wake lean in wake width : 

-5080 do not account for streamtube contraction : 

-1 in cascade parameters : 

-5082 use reynolds stress tensor to compute : 

-1 upwash (rather than TKE) : 

-5084 do not check to see if cfd or experimental : 

-1 data covers one full pitch : 

9999 end of author file 


#1 

#3 


hub 

tip 


1.00 

1.00 

1.00 

0.62 

0.62 


C' 
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Appendix K. — Directory Structure of ~/Codes/fans_nasa 


| ~/Codes/fans nasa/ 



bfans/ 

lib/ 

lib/graf ic/ 

setup bfans/ 

Makefile . Ibm6000 

libgraf ic . a . Ibm6000 

Apollo . m4 

grf lip ,m4 

Makefile . Ibm6000 

Makefile .Iris 

libgraf ic . a . Iris 

Athena . m4 

grgrey . m4 

Makefile .Iris 

Makefile . Sun4 

libgraf ic . a . Sun4 

Aviion .m4 

grgril ,m4 

Makefile . Sun4 

Makefile . wnt 


Dec .m4 

grgrid.m4 

add points. F 

angles . F 

makef lags/ 

Hp9000 . m4 

grinit ,m4 

angles . F 

author . def 

makef lags . Ibm6000 

Ibm6000 ,m4 

grinpa . m4 

arrowhead. F 

bbcascade . F 

makef lags . Iris 

Iris .m4 

grinpf .m4 

author . def 

bf author. m 

makef lags . Sun4 

Makefile . Ibm6000 

grinpi .m4 

cal coeff.F 

bfans ibm6000.F 

makef lags . wnt 

Makefile . Iris 

grkeyl .m4 

calc bl.F 

bfans iris.F 


Makefile . Sun4 

grkey2 . m4 

calc rotor wake.F 

bfans sun4.F 


Newsys ,m4 

grkey3 . m4 

cascade . F 

bfans wnt.F 


Stellar .m4 

grklil .m4 

chan duct geom.F 

capz . F 


Sun4 . m4 

grklin .m4 

circum average. F 

cascade2 . F 


Titan .m4 

grklst ,m4 

constants . def 

constants . def 


Vapor .m4 

grkpla .m4 

converts . def 

converts . def 


Vax .m4 

grlinl ,m4 

cross ref.F 

cross ref.F 


XApollo . c 

grline ,m4 

debug add points. F 

dateandtime . F 


XAthena . c 

grloc2 ,m4 

debug calc bl.F 

error . F 


XAviion . c 

grloct ,m4 

debug calc rotor wake.F 

exact sol hub.F 


XDec . c 

grmesg.m4 

debug estimate. F 

f ex . F 


XHp 9 0 0 0 . c 

grmode . m4 

debug geom.F 

f ex2 . F 


XIbm6000 . c 

grmon2 . m4 

debug get bl.F 

f plus minus2.F 


XIris . c 

grmon3 . m4 

debug get duct geom.F 

fbess j . F 


XNewsys . c 

grmovl . m4 

debug get flow data.F 

fbess j 0 . F 


XStellar . c 

grmov2 . m4 

debug get geom.F 

fbess j 1 . F 


XSun4 . c 

grmov3 . m4 

debug get rotor wake.F 

fbessy.F 


XTitan . c 

grmovy .m4 

debug get row geom.F 

fbessyO . F 


XVapor . c 

grnear . m4 

debug intersect.F 

fbessyl . F 


contng.m4 

grperl ,m4 

debug normalize. F 

fgen . F 


conxch . m4 

grpers . m4 

debug scales. F 

get alpha range. F 


gr2alst .m4 

grplst ,m4 

debug stream line.F 

get freqs.F 


gr2axes ,m4 

grply2 . m4 

disp . F 

get input. F 


gr2ctrl ,m4 

grply3 . m4 

error . F 

get nu range. F 


gr2curs ,m4 

grprnt . m4 

estimate2 . F 

get pratios reba.F 


gr2get ,m4 

grscal .m4 

flow . F 

getalphas . F 


gr21inl ,m4 

grscpt ,m4 

geom.F 

getdeltas . F 


gr21ine ,m4 

grsfac .m4 

get bl.F 

input bb.F 


gr2plst ,m4 

grsset ,m4 

get duct geom.F 

input new.F 


gr2set ,m4 

grsymb .m4 

get file name.F 

interp . F 


gr2sfac ,m4 

grthrl ,m4 

get flow data.F 

jsp te.F 


gr2sset ,m4 

grthre ,m4 

get geom.F 

jsplit.F 


graf ic . doc 

grtime ,m4 

get mesh data.F 

lusolve . F 


gralst ,m4 

grtran .m4 

get rho sound. F 

modes . F 


granot . m4 

grvalu ,m4 

get rotor wake.F 

mysqrt . F 


graxes . m4 

grvecl ,m4 

get row geom.F 

nasa s.F 


grcf il ,m4 

grvect ,m4 

get velocity data.F 

parameters . def 


grcler ,m4 

mrconl .m4 

interp . F 

parameters . fgen . def 


grclos ,m4 

mrcont ,m4 

intersect2 . F 

power ratio. F 


grcnew . inc 


leasqr . F 

read author. F 


grcntr ,m4 


ludcmq. F 

read input. F 


grcol2 ,m4 


mkdir matlab.F 

response . F 


grcol3 ,m4 


naca0012 . F 

run bbcasc.F 


grcolr .m4 


normalize . F 

s pred.F 


grcomn . inc 


parabola . F 

s spec.F 


grconl .m4 


parameters . def 

trapezoid. F 


grcon2 .m4 


plot bl.m 

turbspect . F 


grcont . m4 


plot bv.F 

zeros ex sub.F 


grctrl ,m4 


plot geom.F 



grcube . m4 


plot wake.m 



grcurs . m4 


read author. F 



grcutt ,m4 


reduce .m 



grdash .m4 


save get bl.F 



grdrw2 .m4 


sbf author. m 
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Appendix L. — SetupBFaNS Subroutine Hierarchy 

Bold items are called by multiple routines 


setup_bfans 


1. error 

A. 

read author 

L. 

triangle 


1. 

cross ref 


1 . cal_coeff 

B. 

mkdir matlab 


2. vect_comp 

C. 

get_geom 


3. plot_bv 


1. 

disp 


a. grinit 


2. 

get_duct_geom 


b. arrowhead 



a. getfllename 


c. grscpt 



b. error 


d. grline 



c. sort 

M. 

geom 



d. debug_get_duct_geom 


1. error 


3. 

get_row_geom 


2. parabola 



a. getfllename 


3. debug_geom 



b. error 

N. 

flow 



c. debug_get_row_geom 


1. error 


4. 

debug_get_geom 

0 . 

cascade 

D. 

disp 


1. error 

E. 

get_flow_data 

P. 

naca0012 


1. 

get file name 

Q- 

angles 


2. 

get mesh data 

R. 

upwash 



a. error 

S. 

get_bl 


3. 

get_velocity_data 


1. error 



a. error 


2. interp 



b. stress 


3. save_get_bl 


4. 

circum_average 


4. debug_get_bl 



a. error 

T. 

calc bl 


5. 

get_rho_sound 


1 . debug_calc_bl 



a. error 

U. 

get rotor wake 



b. interp 


1. error 


6. 

stream function 


2. interp 



a. error 


3. debug_get_rotor_wake 


7. 

debug_get_flow_data 

V. 

wake lean 

F. 

add_points 


1. leasqr 


1. 

error 


a. ludcmq 


2. 

debug_add_points 


b. solnq 

G. 

normalize 

w. 

calc rotor wake 


1. 

error 


1. error 


2. 

debug_normalize 


2. interp 

H. 

stream line 


3. debugcalcrotorwake 


1. 

error 

X. 

chan_duct_geom 


2. 

interp 


1. error 


3. 

debug_stream_line 


2. plot_geom 

I. 

intersect2 


a. grinit 


1. 

error 


b. interp 


2. 

interp 


c. error 


3. 

debug_intersect2 


d. gr2sset 

J. 

estimate2 


e. gr21ine 


1. 

error 

Y. 

scales 


2. 

interp 


1 . debug_scales 


3. 

debug_estimate2 

Z. 

write file 

K. 

circumaverage 


1. error 
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Appendix M. — Cascade Parameters 

For a two-dimensional cascade with incompressible flow, the cascade loss, drag and lift coefficients 
coefficient can be determined from, 


co = - 


cosp 


te 


(8 2 \ 

f-ll 

r cos P/ e ^ 

lb f 


v COS Pfe y 


r- _ X cos 3 

C D=™~ 2F“ 

b cos 7 


C L =2^(tanp fe -tanP ?e )cosP„, -C^tanP,, 

b 


P„,=tan -.f tan Pfe^ ta °P« 


(see Ref. 11, Eq. (3.3)2) 
(see Ref. 11, Eq. (3.33)) 
(see Ref. 11, Eq. (3.18)) 
(see Ref. 11, Eq. (3.6)) 


These incompressible cascade relationships assume that the axial velocity at the trailing edge of the 
cascade is the same as that at the leading edge. However, in a real engine, the axial velocity may change 
significantly across the cascade due to contraction in the flow path. The following derivation attempts to 
include the effect of streamtube contraction on the cascade parameters. The derivation assumes steady, 
incompressible flow through a stationary, rectilinear cascade. The flow is assumed to be uniform 
upstream and downstream from the cascade. 


M.l Total Pressure Loss 

The cascade loss coefficient will be defined as, 


co = — 


f i \ ( \ ^ 

Pi + -p Wl 2 - p 2 + -pw? 

z y v 2 


1 


pWi 2 


n 



t i zr r 




i _L i J_M J_ i_ J. 


w v . 


\ 

5, 


w v . 


/ 

r= 

i*- p 

J+. 

i 

i 


t t t t P ,t t t t 


Figure M.l. — Cascade control volume. 
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From the loss coefficient, the pressure difference across the cascade is, 


^-P 2 =lp[(oi-l)lF , 2 + lF 2 2 ] 


PI-P2=~P\ 


(co-i) 


W 


xA 


COS Pi 


f W, - v 


x,2 


V. C 0 S P 2 J 


P t -Pl=\ pw x 2 , 


(s-i(— 

i^cos 


A 2 ( nr \ 2 { 


Pi 


W, 


x,2 


yW Xy i 


v cosp 2 J 


If we define 7 = ff r _2 /W x , 1 , then 


Pl _P 2= 


(co-l) 


— T 

v c°spi y 


7 V 

X 

v COS P2 , 


M,2 Lift and Drag Definition 

We will define the mean velocity and mean angle such that, 


W - 

rr m 


W Xtl 

cosP„, 


The equation used to determine fi m will be derived later. We will define lift as the force component 
that is perpendicular to the mean direction, and drag as the force component that is parallel to the mean 
direction. 



x 

Figure M.2. — Lift and drag forces. 


N ASA/CR— 20 1 0-2 1 6898/VOL 1 


58 



The lift and drag are related to axial and tangential forces via, 


'd' 


r cosP m sin P m ' 

(X' 

,L) 


v - sin Pm COS p ;tt , 

J, 


or inversely 




COS Pm -sin Pm" 


J, 


v Sin p,„ COS pm , 



We will define the lift and drag coefficients as, 

C L = -y — and C D = y — - 

2 PW " bh ' 2 P ’ r " bh ' 

The “hats” on Cl and C/j indicate that the lift and drag are defined as being positive in the x' and / 
directions shown in Figure M.2 In reality, the lift and drag forces on the control surface are in the 
opposite direction to those shown in Figure M.2. Therefore, the actual lift and drag coefficients (without 
the “hats”) are: 


C L = 


-L 


]_ 

2 


P W%bh\ 


and Cl, 


-D 

\pW%bh x 


We will define the axial and tangential force coefficients as, 


Cx 


X 


P w I,C h \ 


and Cy - 


pW^xh 


The lift and drag coefficients are related to the axial and tangential force coefficients via, 


- c D ' 


( W x,l ) 

^ / 

T 

k~ C L 2 


{ w m J 





( ]y 

rr m 

b_( 

SY L 


l^.i J 



cosp m sin(3 w Y Cx 

-sinP,„ cosP m Jy Cy 
or inversely 

cos P W j -sinP w Y-C D 
sin Pm cosp m X~ C L 




J 




J 


M,3 Axial Momentum Equation 

£F x = ffr xP tr-dA 

A 


x+p 1 a 1 -p 2 a 2 -p{a 1 -a 2 )= W x , 2 (p W x 2 A 2 ) - W x \ (p w xJ 4 ) 


From mass conservation, we know thatpfV x 2 A 2 - p W X \A\ . Therefore, 
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X + P l A l -P 2 A 2 -P(4-i 2 )=PMfe-^i) 
x = p 2 a 2 - p,a, +p(a 1 -a 2 )+ p wl j4 (x - 1) 

The average pressure on the endwalls will depend on the endwall geometry. For simplicity, we will 
assume that, 

P=\{Pl+Pl) 


Thus, the x-momentum equation becomes, 


x = ^p 2 -pMi 


( A A 

1+2^ 


A 


i ) 


+P^Vt(x-i) 


After substituting the pressure difference from the loss equation, we get, 


X = --pW x ^A t 


(ffl-l) 


cospj 


+ 


COSp2 


f a A 

i 4 

4; 


+ i 


P^ a 4(x-i) 


Therefore, the axial force coefficient is. 


Cx = 


X 


1 A 


\pWhxA 2 t/7 i 


(co-l(— 

ycos 


' 2 ' r * 


Pi 


+ 


cosp 2 


f A ^ 

1 + A 

4 


+ 2^-(l- X ) 

xh\ 


Since /(| = x/?| , then. 


Cx - 


X 


p^ll/q 


(c° — 1 ) 


A 2 


COS Pi 


+ 


A 2 


COSP2 


f A A 

1 + 2^ 
A 


+ 2(l-x) 


From mass conservation for incompressible flow, the area ratio is related to the velocity ratio via, 

A 2 _ W x> i _ 1 


A w Xt2 x 


Therefore, the axial force coefficient becomes. 


Cx=~- 


(co-lf — — 

\ COS Pi 


A 2 f A 2 

X 


COSp 2 


'1+1 

V X , 


+ 


2(1 -x) 


A 2 


A 2 


Cx =- 


cos Pi 


COSp 2 _ 


'l+x' 

v X 


+ 


2 (l-x)-T 


CO 


f i Y 

T 1+x ] 

(cos Pi J 

l X J 
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M.4 Tangential Momentum Equation 


Y J Fy=§W y pW -dA 

A 

Y = W ya (pW xa A 2 )-W y ,i(pW xA A l ) 

From mass conservation, we know that pW x jAj - pW x \A\ . Therefore, 

Y = pW xJ A l (jV ) , 2 -Wyj) 

Y = P fY x ,iA (\V X>2 tan p 2 - W xA tan p t ) 

Y = P fF x 2 i4 (x tan p 2 - tan pi ) 
Therefore, the tangential force coefficient is, 

Cy =-j = 2(x tan p 2 - tan Pi ) 


M.5 Mean Cascade Angle 

Substitute the expressions for C x and C Y into the equation for C D , 


-c D = 


V tfm j 


1 


\2 


COS Pi 


X 


\2 


cosp 2 


^i + x A 


V ^ J 


+ 2(l-x)-^®| 


\2 


COsPi 


^1 + X A 


V ^ J 


cosp,, 


+ 


2(x tan p 2 - tan Pi )sin P, : 


We will now choose the mean angle such that the drag coefficient is directly proportional to the loss 
coefficient. This requires that the following condition be satisfied, 


1 


f 1 1 

2 

f X | 

2" 

(i + x) 

+ 2(l-x) 

2 


l cospi. 


lcosp 2 . 


{ x J 


cos 


P„, + 2(x tan p 2 - tan Pi )sin p,„ = 0 


This can be solved for P„,to give the following expression, 


tan P,„ = - 


cosPi 


cosP 2 


ri±x^ 

v X 


+ 


(i-x) 


tanp t -%tanp 2 

After a little algebra, we can show that for x = 1 , the above expression reduces to the equation shown 
earlier for a 2-D cascade, 

tan P„, = i (tan p 2 + tan Pi ) 
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M.6 Drag Coefficient 

The condition above also results in the following expression, 


r 1 — 1 
-Cn = — co — 

2 b 


1 


COS Pj 


(t. nr ^ 2 


cosP„ 


W-, 


xA 


V ^’i)i J 


\ x 


w xl 

Since W m = — - — , then, 
cosp m 


C D = - co- 


1 _ T COS 3 P 


2 b cos 2 Pi 


1 + X 
l X 


For x = 1, the above expression reduces to the equation shown earlier for a 2-D cascade, 

rccs’A, 


b cos J3 l 


M.7 Lift Coefficient 


Cy 




\W X , l 


- (- C D sin p m - C L cos p,„ ) 

T 


-c, =- 


cosp,, 


fvtr ^ 2 


w, 


xA 


V J 


— Cy + C [) Sill [ f 
b 


W . i 

Since W m = — — — , then, 
cosP„, 


~C L = cos P„, \C Y +C D tan p,„ 
b 

~ C L = 2 y (x tan p 2 - tan Pi )cos P,„ + C D tan p,„ 
b 

C L =2^(tanp 1 -x tan P2)cosP,„ -C c tanP m 
b 

which is the same as that for a 2-D cascade when x = 1. As noted in Reference 1 1 (p. 59), the drag 
contribution to the lift coefficient is very small; therefore, the C D tan P„, term is neglected in 
SetupBFaNS. 
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Appendix N. — Reynolds-Stress Transformation 

N.l Generalized Tensor Transformation 

Any coordinate transformation can be written in the following form (Ref. 12, p. 31), 


dvi = dx ■ 

' dx j 

where Xj is the original coordinate system, and y, is the new coordinate system. If we limit ourselves to 
rectangular Cartesian coordinate system (which can only undergo rotations, reflections, and translations), 
then each term in the transformation matrix is constant, and we can integrate to get, 

yi = P yXj + a. 

The coefficients, (3,y, are the direction cosines between the new coordinate axes, y„ and the old coordinate 
axes, Xj. 

For any second-order tensor (such as Reynolds stress), the transfoimation from one coordinate system 
to another can be written as (Ref. 12, p. 33), 


rr' _ -r fyi 8y j 
1 ij ~ 1 mn _ - 

OX m OX n 

For rectangular Cartesian coordinates, this transformation becomes, 

Tij = T’mnP/mP. jn 

where Tj is the tensor in original coordinate system, and T{j is the tensor in the new coordinate system. 


N.2 Transformation from Engine Coordinates to Duct Coordinates 

Flanson derived his turbulent inflow noise theory for a rectilinear cascade (Ref. 4). Fie defined the 
duct coordinates such that x d is in the axial direction, y d is in the tangential direction (y d = r9), and z d is 
binoimal to x d and y d . Fie assumed that the velocity was uniform and that the binormal component was 
zero. 

To satisfy the assumption that the binormal component is zero, Setup BFaNS defines the duct 
coordinates such that x d is aligned with the meridional flow direction, y d is in the tangential direction 
(y d = r9), and z d is binoimal to x d and y d . If we denote the duct coordinates by (x d , y d , z d ) and the engine 
coordinates by (x e , y e , z e ), then the transfoimation from engine coordinates to duct coordinates is, 




r cos(90-<))) 

0 

sin(90-<))) N 


yd 

= 

0 

1 

0 


\ Z d V 


v — sin(90 — <()) 

0 

cos(90-<()) y 

v Z <?y 


where ()> is the meridional flow angle. Therefore, the transformation matrix simply becomes, 



' sin 4> 

0 

cos 



0 

1 

0 



v — COS ()) 

0 

sin 
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N.3 Reynolds Stress-Transformation 

As noted earlier, the transformation for a second-order tensor is, 


We will let Ty represent the Reynolds stress (c,c ,• ). Since the Reynolds-stress tensor is symmetric, so we 
only need to consider the following six components, 

T\\ = T\ 1P1 ipi i + T22P12P12 + T33P13P13 + 

+ -^1 3 (P 1 3 P 1 1 + Pi lPl3 ) T23 (P 1 3 P 1 2 +P12P13) 

2 ~ T \ 1P21P21 +T22P22P22 +T33P23P23 + T12 (P22P21 +P21P22) 

+ -^1 3 (P23P2 1 +P21p23) +? 23(P23P22 +P22P23) 

T33 = T \ 1P31P31 +T22P32P32 +T33P33P33 + ^12 (P32P3 1 +P31P32) 

-^13 (P33P3 1 + P3 lP33 ) + ^23 (P33P32 +P32P33) 

T\2 = W11P11P21 +T22P12P22 +T33P13P23 + T\2 (P12P2I +P11P22) 

+ -^1 3 (P 1 3 P 2 1 + Pi 1P23 ) ^ 23 (P 1 3 P 22 +P12P23) 

T { 3 =T 11 p 11 p 31 +722P12P32 +T33P13P33 +7i 2 (Pi2P3i +P11P32) 

+ -^1 3 (P 1 3 P 3 1 + Pi lP33 )+ -^23 (Pi 3P32 +P12P33) 

^23 = T11P21P3I +T22P22P32 +T33P23P33 + -^1 2 (P 22 P 3 1 +P21P32) 

+ T| 3 (P 23 P j \ +P2lP33) +? 23(P23P32 +P22P33) 

For the conversion from (x e , y e , z e ) to (x d , y d , z d ), we get, 

7 |'| = T\\ sin 2 ()) + 733 cos 2 ()) + 2^3 sin <|) cos (p 
T22 = T 2 2 

J33 = T\ 1 cos 2 (p + 733 sin 2 (p — 27 j 3 sin <|) cos (|) 

T{ 2 = T\2 sin i|) + 723 cos 

71 ' 3 = - 7 j 1 sin <p cos (|) + 733 sin (p cos (p + 7 ] 3 (- cos 2 (p -1- sin 2 cp) 

?23 = - 7|2 cos 4 + ?23 sin (|) 
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Appendix O. — Upwash Computation 


In Reference 4, Hanson derived the transformation that converts from duct coordinates to cascade 
coordinates, 



( 

f n 


( K ' 


f n jT 


^ n ^ 

\ 



cos| 

--A 

cosn/5- + sin 

hr -5 

Jsinvj/^ sinvj/^ sin 
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sin\|/^ sinij/^ -cos\\r L sini}/^ 





U J 

/ 

v 2 j 


1.2 j 

/ s 
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. n y 
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\ 
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-sin c 

l 2 

; cosy^ 


cos 

cos\j/£ -sin\|/^ 

yd 

\ Z C J 



(n , A 


r n ^ 


1 71 

\ ^ 

1 • [ 

' n S\ 


\ Z d j 


cos 

•i 

1 

JA 

sin\|/^ -sinj 


sin\|/^ cosij/^ sin 

'I 

1 

siny 5 +cos 

•i 

1 

sin\|/^ cos\|/5 cos\j/^ cos\j/^ 



V 



V 2 J 



\ 

V J 

y 



where x d is in the axial direction, y d is in the tangential direction (y d = r9), z d is binormal to x d and y d . Due 
to the difference in stagger conventions between Hanson’s derivation and SetupBFaNS, Hanson’s 
stagger angle has been replaced with the complement of q. 

Hanson’s derivation assumed that the binormal velocity component is zero. Therefore, in Hanson’s 
derivation, the geometric lean and sweep are the same as the aerodynamic lean and sweep. However, 
Setup BFaNS defines the duct coordinates such that x d is aligned with the meridional flow direction 
rather than the axial direction. Thus, in Setup BFaNS, the geometric lean and sweep are not the same as 
the aerodynamic lean and sweep. To account for this difference, we must use aerodynamic sweep and 
lean in the transfoimation matrix, i.e., 




|cos\|/5 + sin|^— - £, jsiny^ siny's 

-sinj^-^Jcosyi 

k |sin\| i's -sinf-^-Msinyi eosy^ 



cosy 5 - cos 



sin y 5 + cos 


^-^jsiny£sinys 

cosy), 

^-Usinyi, cosy's 


\ 

-cosy) sin y s 
-siny) 
cosy) cosy) 


The derivation for aerodynamic sweep and lean is shown in Appendix P. The upwash is defined as the 
Reynolds-stress component that is normal to the airfoil surface (i.e., in the y c direction). Therefore, the 
upwash can be determined from, 


T7.2 ~ Tmnfilmfiln 


where T’ mn is the Reynolds-stress tensor in the duct coordinate system (i.e., x d , y d , z d ). 
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Appendix P. — Aerodynamic Sweep and Lean 
P.l Coordinate Transformation 

Hanson derived his turbulent inflow noise theory for a rectilinear cascade (Ref. 4). He defined the 
duct coordinates such that x d is in the axial direction, y d is in the tangential direction (y d = r9), and z d is 
binormal to x d and y d . He assumed that the velocity was uniform and that the binormal component was 
zero. 

To satisfy the assumption that the binormal component is zero, Setup BFaNS defines the duct 
coordinates such that x d is aligned with the meridional flow direction, y d is in the tangential direction 
(y d = r9), and z d is binormal to x d and y d . If we denote the duct coordinates by (x d , y d , z d ) and the engine 
coordinates by ( x e , y e , z e ), then the transformation from engine coordinates to duct coordinates is, 

f cos(90-<))) 0 sin(90-<t>)Y x e ' 

yd = 0 1 0 Te 

K z dJ y— sin(90 - <)>) 0 cos(90 - z e y 

where <j> is the meridional flow angle measured from the z d axis. Therefore, the transformation matrix 
simply becomes, 

^ sin (() 0 cost)) 

h= 0 10 

y -cos(|) 0 sin 4> 

Thus, the relationship between duct coordinates and engine coordinates is, 

Xd = x e sin ( j> + z e cos ()> 

yd =y e 

Zd = —x e cos <j> + z e sin <)> 

P.2 Geometric Sweep and Lean 

The engine coordinates of the airfoil leading edge are defined parametrically as a function of z e , 

x e,le ~ ^e,le\ z e) 
ye,Ie ~ Te,/e( z e) 

From this parametric representation, we can determine the geometric sweep and lean via, 

dx e j e 

tanys =— — 

a^e.le 

. dy e j e 

tan \| ) L = —^— 
d z e Je 

P.3 Aerodynamic Sweep and Lean 

Aerodynamic sweep and aerodynamic lean are defined in the (x d , y d , z d ) coordinate system as, 
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, r d x d,le 

tan \\i s = 

dz d,le 


, we see that, 


From the coordinate transformation, 

dx e j e . dz e i e 
tan\|/£ = smqH — cosq) 

d z d,le dZgjjg 


After applying the chain rule, we get, 


, \ dy e j e 

tan\) f L = —^— 

dZd,le 


dxp i p . 

— smi|) + cos(|) 

v dz e i e 


i uz e i e J 

tan\|/ 5 = -A— 

ClXp Ip i . i 

— cos (j) + sin (j) 

dzgig 


(a \ 

“Ve.fe 

v dz e ]g 


tan \| n 

ClXp Ip i . i 

— cos([) + sm(t) 

v dz e i e y 

After substituting the definition of geometric sweep and lean, we get, 

, tan w c sin cb + cos cb 

tan \\f' s = — — 

-tani^ cos t)> + sin tj> 

tan\\) L 


After a little algebra, we finally get, 


, tan vi i j 

tanv|/£ = — 

-tanv|/£ cos c)) + sin tj> 


Vs =90 + \|/ s -c|) 
tani|/£ 


, tan vi/ j 

tan\| f L = — 

- tan v^ cos tj> + sin c() 
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Appendix Q. — Boundary-Layer Estimation 

Q.l Manually Specifying Boundary-Layer Thicknesses 

The boundary-layer thicknesses at Stations 1 and 3 can be manually specified in the SetupBFaNS 
input file. If the boundary layers are not manually specified, one of the following procedures will be used 
to estimate them. 


Q.2 When the Flow Field is Determined from an Inviscid Flow-Field Prediction 

For an inviscid flow-field prediction, Setup BFaNS assumes a flat-plate turbulent boundary layer. 

The boundary-layer thickness is computed from, 

p Ua,L 

‘<L = 

V 

5 = 0.37fr^ 1/5 (Ref. 13, p. 638) 


where U x is the free-stream velocity relative to the surface, and L is the distance over which the boundary 
layer grows. The kinematic viscosity is hardcoded to 1.5x10° m 2 /s (air at 20 °C — see Ref. 8). The free- 
stream velocity is estimated from the flow speed at the surface, 


U a 


= r 2 + c 2 

y ' l '~r,surf 


+ 




where 


U S urf = rotational speed of the surface 
The distance over which the boundary layer grows is estimated from, 

L = r Ux ( z - z o) 

'~'z,surf 


where 


z 0 = origin of the boundary layer 

The following table shows the parameters used for z 0 at each station. 


Location 

Value used for Zq 

Station 1, Hub 

7 

nose 

Station 1, Tip 

Z hish 

Station 3, Hub 

Z core 

Station 3, Tip 

Zhish 


Note that the boundary-layer origin for “Station 3, Tip” is assumed to be the hilite. This can be a gross 
approximation since the flow has passed through the fan. 

Based on turbulence profiles shown in Reference 14, the average turbulence in a flat-plate turbulent 
boundary layer is approximately 5 percent of the free-stream velocity. For simplicity, SetupBFaNS 
assumes a uniform 5 percent turbulence intensity across the boundary layer. 
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Q.3 When the Flow Field is Determined from a Viscous CFD Prediction 

In classical boundary-layer theory, the boundary-layer thickness is defined as the point where the 
velocity is 99 percent of the free-stream velocity (Ref. 8, p. 388). Unfortunately, the velocity profiles in 
the fan stage are very non-uniform, making it difficult to determine a “free-stream” velocity. 

A method was developed to estimate the boundary-layer thickness based on the circumferentially 
averaged turbulent kinetic energy profiles. Away from the endwalls, these profiles are much flatter than 
the velocity profiles, making them a better discriminator for determining the edge of the boundary layer. 
An example of this procedure is shown in Figure Q. 1 . 

At each axial station, SetupBFaNS divides the flow field into a “tip region” and a “hub region”. The 
radial location that corresponds to the minimum turbulent kinetic energy divides these two regions. 

Setup BFaNS computes the slope of the TKE profile in each region, and normalizes the slope such that 
its maximum value is one. Setup BFaNS then searches for the point where the normalized slope is equal 
to a tolerance value, which is hardcoded to 0.1. This point is defined as the edge of the boundary layer. 
Via interpolation, Setup BFaNS then computes the free-stream velocity (relative to the surface), the free- 
stream meridional velocity, and the turbulence level at the edge of the boundary layer. 


Station 3 - Hub Station 3 - Tip 






Figure Q.l. — Using K to determine edge of boundary layer. 
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Q.4 When the Flow Field is Determined from Experimental Data 

For experimental data, SetupBFaNS uses the same procedure that it uses for a viscous CFD 
prediction. 

Q.5 Extrapolation to Blade and Vane Leading Edges 

Once the boundary layers are determined (or manually specified) at Stations 1 and 3, they are 
extrapolated to the blade and vane leading edges. The extrapolation assumes a flat-plate turbulent 
boundary layer. 
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Appendix R. — Rotor-Wake Estimation 

R.l When the Flow Field is Determined From an Inviscid Flow-Field Prediction 

For an inviscid flow-field prediction, SetupBFaNS uses empirical correlations to estimate the rotor 
wake properties. In Reference 15, Glegg used the following correlations to estimate the wake velocity 
defect (AW), wake thickness (5,,,) and wake turbulence (w'w'), 


A W = - 


K l 

\w\ 

2 

s 2 

s - 

-Ki8 2 


K = 2^82(5-^62) 


w'w' = K 3 (AW) 2 

These correlations are based on work by Wygnanski et al (Ref. 16). The constants depend on the shape of 
the body that generates the wake. Following Glegg’s example, Setup BFaNS uses the constants that 
correspond to a symmetric airfoil, which are, 


K 0 =1.56 
K\ = -380 
K 2 = 0.32 
K 3 =0.15 

Based on these correlations, the only information that is needed to compute the wake properties is the 
momentum thickness, the tree-stream velocity, and the streamwise distance downstream from the airfoil 
trailing edge. 

As shown in Section 4.9.4, the momentum thickness can be determined from, 

^ = 0.006 + 0.0002 e 7 - 5D f _1 

b 

If the axial velocity changes between the blade trailing edge and Station 3, then the momentum thickness 
is adjusted according to, 


Station 3 


r w ZJe ^ 2 

(_ W z,3 


The tree-stream velocity is assumed to be the circumferentially averaged relative flow speed at Station 3. 
The streamwise spacing between the fan trailing edge and Station 3 is, 


W - 3 / x 

5 = (z 3 -Zfe) 


W- 


z, 3 


R.2 When the Flow Field is Determined From a Viscous CFD Prediction 

For viscous-CFD predictions, Setup BFaNS determines the 0-locations where, 
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W-W 

vv vv min 
Rmax — ^min 


0.5 


These locations define the edge of the wake. The wake thickness is defined as, 

5 w = r(e+-e-) 


Figure R.l demonstrates this process. 

[ 2 

The wake turbulence is defined as the maximum value of J — K that occurs in the blade passage. It is 
assumed to be uniform across the width of the wake. 

R.3 When the Flow Field is Determined from Experimental Data 

For experimental data, SetupBFaNS uses the same procedure that it uses for a viscous CFD 
prediction. 

R.4 When the Flow Field is Determined from Experimental Data 

For experimental data, Setup BFaNS uses the same procedure that it uses for a viscous CFD 
prediction. 



y/ T 

Figure R.l . — Wake terminology. 
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Appendix S. — Estimating Integral Length Scale 

S.l Integral Length Scale From Empirical Correlations 

Endwall Boundary Layer . — In Reference 15, Glegg used a Von Karman type spectrum to model the 
boundary-layer turbulence interacting with the fan tip. In his model, Glegg used the following 
relationship, 


**8 = 1.2 

where the turbulence wavenumber (k e ) is related to the integral length scale via, 

_ r(5/6) Vtt _ 0.747 
e ” r(l/3) A _ A 

After substituting Equation (S2) into Equation (SI), we get the following relationship, 

A = 0.625 


(SI) 


(S2) 


Following Glegg’ s example, we will assume that the integral length scale is constant across the layer, and 
is proportional to the thickness of the layer, 

A = C5 


The default value for the coefficient (C) is 0.62, consistent with the relationship that Glegg used. The user 
can over-ride this default with Card 5062 (Station no. 1, ID), Card 5064 (Station no. 1, OD), Card 5066 
(Station no. 3, ID) and Card 5068 (Station no. 3, OD) in the author file. 

Rotor Wake . — In Reference 15, Glegg also used a Von Karman type spectrum to model the wake 
turbulence interacting with the FEGV. In his model, Glegg used the following relationship, 

k e L 0 =1.81 (S3) 


where 



After substituting Equation (S2) into Equation (S3), we get the following relationship, 

A = 0.685 w 


Following Glegg’ s example, we will assume that the integral length scale is constant across the wake, and 
is proportional to the thickness of the wake, 


A = C5 W 


The default value for the coefficient (C) is 0.68, consistent with the relationship that Glegg used. The user 
can over-ride this default with Card 5070 in the author file. 

S.2 Integral Length Scale From Turbulence Model 

The integral length scale is assumed to be proportional to the turbulent mixing length, 
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A = C 


3 

K2 


8 

The default value for the coefficient (C) is one. However, the user can over-ride this default with Card 
5056 in the author file. At the walls, the turbulent length scale is hardcoded to be zero. 


S.3 Circumferential Averaging 

Results from Reference 17 showed that the Liepmann spectrum provided a good fit to the 1-D 
turbulence spectra measured downstream from a 22 in. fan rig. In addition, the results showed that the 
circumferentially averaged axial, tangential and radial turbulence components were nearly equal. Those 
observations suggest that the turbulence may be modeled with an isotropic spectrum model, even though 
the turbulence is non-homogenous. 

We will assume that the Liepmann spectrum adequately represents the circumferentially averaged 
turbulence. Therefore, the 3-D upwash spectrum is (Ref. 4), 


4*22 y 


Ik. 

V- 


avg 


A5 


avg 


1 + A ~ aV g 


k r + A; 2 


(A j- + A 2 -t- 


In this equation, K avg and A avg represent the turbulent kinetic energy and integral length scale associated 
with the circumferentially averaged turbulence. 

Reference 1 7 describes a procedure for determining K avg and A avg from experimental data. When 
experimental measurements of the integral length scale are not available, Setup BFaNS uses the 
procedure described below. 

We will assume that the turbulence spectrum at each point can be represented by the Liepmann 
spectrum. With that assumption, the longitudinal turbulence spectrum at each point is (Ref. 18), 
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where K and A depend on (r, 0, z), and k\ = 00 / C . If we circumferentially average the longitudinal 


spectrum, we get, 
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At kx = 0, the average longitudinal spectrum evaluates to, 
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Since we have assumed that the average spectrum can be adequately modeled with the Liepmann 
spectrum, then 
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At k\ = 0, this becomes, 
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After equating (S4) to (S5), we find that, 
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Thus, the average length scale is, 
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